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PREFACE 


The  report  herein  was  sponsored  by  Headquarters,  US  Army  Corps  of  Engi¬ 
neers  (USAGE)  and  US  Air  Force  under  the  effort  "Design  for  Heavy-Weight  Air¬ 
craft."  The  USAGE  Technical  Monitor  was  Mr.  Paige  Johnson. 

The  computer  program  was  developed  at  the  US  Army  Engineer  Waterways 
Experiment  Station  (WES)  from  May  1991  through  September  1991  by  the  Pavement 
System  Division  (PSD)  of  the  Geotechnical  Laboratory  (GL) .  This  report  was 
written  by  Mr.  Carlos  R.  Gonzalez,  PSD,  under  the  supervision  of  Dr.  George  M. 
Hammitt  II,  Chief,  PSD;  and  Mr.  Jim  Hall,  Chief,  System  Analysis  Branch.  Tha 
work  was  performed  under  the  general  supervision  of  Dr.  W.  F.  Marcuson  III, 
Director,  GL. 

At  the  time  of  publication  of  this  report.  Director  of  WES  was 
Dr.  Robert  W.  Whalin.  Commander  and  Deputy  Director  was  COL  Leonard  G. 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non- SI  units  of  measurement  used  in  this  report  can  be  converted  to 
SI  (metric)  units  as  follows; 


_ Multiply _ 

Inches 

pounds  (force) 

pounds  (force)  per  square  inch 
square  inches 


-  By 

2.54 

4.448222 

6.894757 

6.4516 


To  Obtain _ 

centimetres 
newtons 
kilopascals 
square  centimetres 
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A  COMPUTER  PROGRAM  FOR  THE  CALCULATION  OF  THE 
EQUIVALENT  SINGLE-WHEEL  LOAD  FACTOR 

PART  I :  INTRODUCTION 

1.  For  many  years  the  Corps  of  Engineers  (CE)  has  been  implementing  the 
concept  of  an  equivalent  single -wheel  load  (ESWL)  to  represent  the  effects  of 
multiple  wheels  in  a  vehicle  gear  assembly  upon  a  pavement  system.  In  the 
early  years  of  the  development  of  the  CE  pavement  design  procedure,  the  Bous- 
sinesq  solution  for  the  effects  of  circular  loads  in  a  linear  elastic  half 
space  was  selected  to  approximate  the  r*.  sponse  of  a  pavement  structure  under 
wheel  loads.  Closed  form  solutions  for  computing  the  pavement  response 
(stresses  or  strains)  under  multiple  wheels  were  difficult  and  tedious  to 
compute.  Even  today  these  solutions  required  considerable  computer  power. 

For  this  reason,  the  concept  of  converting  a  multi-wheel  gear  load  to  a 
single -wheel  load  that  would  produce  the  same  effect  at  a  given  depth  in  a 
pavement  system  was  introduced.  This  conversion  or  equivalency  was  based  on 
vertical  deflection  at  a  given  depth.  Yoder  and  Witzak  (1975)*  and  the  BOEING 
Company  (1963)**  describe  in  detail  the  concept  behind  this  equivalency.  For 
the  purpose  of  this  report,  it  is  sufficient  to  know  that  the  ESWL  is  the  load 
on  a  single  wheel  that  will  produce  the  same  maximum  vertical  deflection  at  a 
specific  depth  as  all  the  tires  in  a  gear  assembly  and  that  the  ESWL  varies 
with  depth.  Figure  1  illustrates  the  ESWL  concept.  The  equivalent  single¬ 
wheel  load  factor  (ESWLF)  is  the  ESWL  expressed  as  a  fraction  of  the  gear 
load.  The  ESWLF  would  be  approximately  equal  to  the  fraction  of  one  tire  load 
t.c  the  gear  load  (tire  load/gear  load)  at  the  surface,  and  approaches  asymp¬ 
totically  an  ESWLF  value  of  1.0  with  -de'pth. 

2.  This  report 'desor-ihes  the  implementation  of  a  computer  program  to 
calculate  the  ESWLF  of  any  gear  configuration  using  Boussinesq  solution  and 
vertical  deflection  as  the  equivalency  parameter.  It  is  not  intended  to  dis¬ 
cuss  the  pros  and  cons  of  the  concept  assumptions,  but  rather  to  present  and 


*  E.  J.  Yoder  and  M.  W.  Witzak.  1975.  "Principles  o&  Pavement  Design," 
Wiley- Interscience ,  NY. 

**  Boeing  Company,  Inc.  1963.  "A  Computer  Program  for  Determining  the  Flo¬ 
tation  Capability  of  Aircraft  Landing  Gear  Using  the  CBP  Method,"  Document 
No.  D64088TN,  Renton,  WA. 
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P1 


P2 


a.  Deflection  under  multiple  gear 


ESWL 


Figure  1.  ESWL  analysis  for  equal  deflection  criteria 
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describe  the  computer  program  developed.  It  has  been  found  that  better  solu¬ 
tions  can  be  achieved  if  more  sophisticated  layered  elastic  system  computer 
programs  are  used,  particularly  when  a  large  number  of  wheels  are  needed  to 
support  heavy  aircraft  or  vehicles.  However,  for  simple  gear  configurations 
with  few  wheels  (i.e.  twin  or  twin-tandem  gears),  the  ESWLF  program  produces 
reasonable  results  and  is  still  in  common  use  even  for  more  chan  four  wheels. 


PART  II:  PROGRAMMING  PHILOSOPHY  AND  SYSTEM  REQUIREMENTS 


3.  Emphasis  on  ease  of  use,  data  entry  flexibility,  and  simplicity  were 
major  concf.ts:iS  while  develc;.ing  this  computer  program.  Menus,  form- like  data 
entry,  and  windows  enhau'  ^  and  simplify  the  manipulation,  calculation,  and 
presentation  of  the  results.  The  computer  code  (Appendix  A)  was  written  with 
the  structure  programming  concept  in  mind  to  improve  code  maintainability. 

The  program  was  compiled  using  the  Microsoft  FORTRAN  5.0  compiler  and  was 
design c;r  to  operate  on  an  IBM  Personal  Computer  or  compatible  machine  under 
the  DOS  3.0  operating  system  and  later  versions. 

4.  The  graphics  libraries  and  extensions  to  the  standard  FORTRAN 


included  with  this  compiler  were  implemented  in  this  program.  Although  not 
required,  an  IBM  AT  class  computer  or  bett^'  with  a  math  coprocessor  and 
640  KBytes  of  main  memory  is  recommended,* 
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PART  III:  USER  INSTRUCTIONS 


5.  The  conunands  and  features  available  in  this  program  are  described. 

An  example  problem  is  used  to  work  the  user  through  the  steps  required  to 
calculate  the  ESWL.  Several  conventions  are  used: 

a.  A  word  or  phrase  inside  brackets  means  that  the  user  needs  to 
type  the  word  or  phrase  at  the  keyboard  or  press  the  indicated 
key.  For  example:  [ENTER]  is  telling  the  user  to  press  the 
enter  key,  [UP  ARROW]  means  pressing  the  up-arrow  key  on  the 
keyboard,  and  so  on. 

b.  The  example  assumes  the  program  is  being  run  from  a  hard  disk 
with  a  logical  drive  C.  If  the  user  is  running  from  a  hard  disk 
other  than  C  or  from  a  floppy  drive,  just  replace  C  in  the 
instructions  with  the  corresponding  drive  letter.  For  the  pur¬ 
pose  of  this  example  the  drive  letter  is  identified  as  C:>. 

This  represents  the  DOS  prompt  and  does  not  need  to  be  typed. 

c.  The  program  supports  the  use  of  the  keyboard  arrow  keys.  In  the 
menus,  the  options  can  be  selected  either  by  moving  a  high¬ 
lighted  bar  to  the  option  desired  or  pressing  the  number  in 
brackets  preceding  the  option. 

Running  the  Program 

6.  Before  running  the  program,  the  hard  disk  or  floppy  disk  where  it  is 
to  be  run  should  be  checked  to  be  sure  sufficient  space  is  available  to  hold 
the  input  and  output  data  files.  To  run  the  program  just  type  C:>  ESWLF 
[ENTER]  at  the  DOS  prompt. 

7.  A  welcome  screen  appears  next  displaying  the  program  title,  author, 
and  date  released.  After  pressing  any  key,  a  menu  containing  the  program 
features  or  options  is  displayed  (Figure  2).  From  this  menu  the  entry  of 
data,  file  manipulation,  and  viewing  of  the  results  can  be  accomplished.  The 
menu  shown  in  Figure  2  allows  the  user  to  build  a  new  data  file,  load  an 
existing  data  file,  save  currently  calculated  results  to  disk,  view  an  exist¬ 
ing  output  data  file,  and  exit  the  program. 


Entering  Data 

8.  In  general,  the  steps  required  to  calculate  the  ESWL  consist  of: 

a.  Gathering  the  vehicle  tire  coordinates,  tire  pressure,  and 
radius  of  tire  contact  area. 
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^  E<iuiwaJ.ent  single  Wheel  L<oad  ^actorCtdiWLit'P  i'rograni  f 


_ ESWLP  PROGHftH  MAIN  HEHU _ 

Cll  Build  a  new  data  file 

C2]  Load  an  existinsr  data  file  from  disk 

C3]  Save  curr«nt  output  data  file  to  disk 

C4]  View  an  output  data  file 

C5]  Exit  pt*ogi*aM 


Figure  2.  ESWLF  program's  main  menu 

b.  Building  the  input  data  file. 

c.  Selecting  the  option  to  proceed  with  the  calculations. 

d.  Viewing  and/or  printing  the  results  obtained. 

9.  The  first  step  in  calculating  the  ESWL  is  to  build  the  input  data 
file.  To  do  so  the  user  would  select  option  number  one  (1)  and  press  [ENTER]. 
The  program  will  then  open  a  window  asking  for  an  input  data  file  name.  For 
simplicity,  the  file  name  extensions  for  the  input  data  file  and  output  data 
file  have  been  preset  to  *.IN  and  *.OUT,  respectively.  At  this  point  the 
desired  file  name,  up  to  eight  characters  long,  can  be  typed  (Figure  3).  If 
the  file  name  already  exists,  the  program  will  warn  the  user  of  such  situation 
and  will  allow  the  user  to  re-enter  file  name.  Pressing  the  [ESC]  key  cancels 
the  input  and  returns  the  user  to  the  menu. 

10.  Once  the  input  file  name  is  entered  a  data  entry  screen  with  the 
required  parameters  that  are  needed  to  perform  the  calculations  is  displayed 
(Figure  4).  A  description  of  the  parameters  in  this  screen  follows: 

a.  Same  contact  pressure  for  all  tires?  (Y/N) 

This  item  permits  the  analysis  of  tires  with  different  contact 
pressures.  However,  it  is  common  for  the  tires  in  a  gear 
assembly  to  be  of  the  same  size  and  have  the  same  inflation 
pressure.  Also,  in  pavement  analysis  the  tire  inflation  pres¬ 
sure  is  often  assumed  to  be  equal  to  the  tire  contact  pressure 
when  the  tire  contact  pressure  is  unknown. 
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Figure  3.  Screen  for  establishing  a  new  input  data  file  name 


Figure  h.  Data  entry  screen  for  the  ESWLF  program 


b.  Same  radius  of  conCacC  area  for  all  tires?  (Y/N) 

This  item  permits  the  analysis  of  tires  with  different  radius 
of  contact  area.  Like  the  case  of  the  tire  pressures,  in  most 
situations  the  radius  of  contact  area  in  a  gear  assembly  is 
also  considered  to  be  equal  for  all  tires. 

c.  Number  of  tires  in  gear  (100  Max) 

This  is  the  nuunber  of  tires  to  be  analyzed.  The  program  allows 
up  to  100  tires . 

d.  Tire  contact  pressure  (psi)* 

If  it  is  decided  to  have  all  tires  with  the  same  contact  pres¬ 
sure,  then  this  item  will  be  active  and  the  value  entered  will 
correspond  to  the  tire  pressure  of  all  tires.  The  tire  pres¬ 
sure  can  be  computed  by  dividing  the  tire  load  by  the  contact 
area  or  it  may  be  assvimed  to  be  equal  to  the  tire  inflation 
pressure . 

e.  Tire  radius  of  contact  area  (in.) 

If  it  is  decided  to  have  all  tires  with  the  same  radius  of 
contact  area,  then  this  item  will  be  active  and  the  value 
entered  will  correspond  to  the  radius  of  all  the  tires.  Since 
the  Boussinesq  theory  only  considers  uniform  loads  on  circular 
areas,  the  radius  of  contact  area  can  be  calculated  by  assuming 
the  contact  area  to  be  circular  and  taking  the  equation  for  the 
area  of  a  circle  and  solving  it  for  the  radius . 

f.  ESWL  radius  of  contact  area  (in.) 

For  the  CE  procedure  this  value  is  always  set  equal  to  the 
radius  of  one  of  the  tires  in  the  landing  gear. 

g.  Number  of  computational  depths  (100  max) 

Because  the  vertical  deflection  caused  by  a  tire  load  varies 
with  depth,  the  program  allows  the  user  to  compute  the  ESWLF 
for  up  to  100  depths . 

h.  Depth  Increment  (in.) 

The  program  uses  the  number  of  depths  from  the  previous  item 
and  the  depth  increment  to  establish  the  depths  at  which  the 
ESWLF  are  to  be  computed. 

i.  Grid  x- increment  (in.)  and  grid  y- increment  (in.) 

Since  this  program  was  designed  to  handle  multiwheel  gear  con¬ 
figurations,  the  maximum  deflection  of  such  a  gear  does  not 
necessarily  occur  under  one  of  the  tires.  Consequently, 
several  points  around  the  gear  must  be  checked  to  obtain  the 
maximum  value.  The  program  does  so  by  building  a  grid  based  on 
the  grid  x  and  y  increments  entered  by  the  user.  The  finer  the 
grid  generated,  the  more  accurate  the  results  will  be.  How 


*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  3. 
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ever,  using  fine  grids  also  increases  the  time  of  the  calcula¬ 
tions.  The  origin  of  the  grid,  which  is  independent  of  the  X-Y 
coordinate  system,  is  automatically  set  to  the  tire  coordinate 
of  the  leftbottom  corner  of  the  gear  assembly.  Figure  5 
illustrates  the  grid  system  the  program  would  generate  for  the 
gear  configuration  shown. 

2-  Symmetry  along  X  and  Y  axes  (Y/N) 

This  feature  takes  advantage  of  gear  configurations  with  symme¬ 
try  along  their  axes.  For  example,  a  gear  is  symmetric  along 
the  X-axis  if  the  tires  on  the  positive  Y-axis  side  and  the 
negative  Y-axis  side  are  mirror  images  of  each  other.  The  same 
analogy  applies  to  the  Y-axis.  This  is  important  because  sym¬ 
metry  greatly  reduces  the  calculation  time.  Figure  5  illus¬ 
trates  the  grid  systems  the  program  would  generate  if  symmetry 
is  not  assumed  along  the  X  and  Y  axes  (Figure  5a)  and  if  symme¬ 
try  is  assumed  along  both  the  X  and  Y  axes  (Figure  5b) . 

Although  the  location  of  the  axes  is  fixed  at  the  center  of  the 
gear,  it  is  allowed  to  set  its  location  at  any  point. 

k.  Enter  or  change  tire  coordinates?  (Y/N) 

This  entry  allows  the  user  to  enter  or  change  tire  coordinates , 
When  answered  affirmative,  a  window  like  the  one  shown  in  Fig¬ 
ure  6  will  appear  to  ask  for  these  values  for  each  of  the 
tires.  The  tire  pressures  and  radius  of  contact  area  are 
entered  only  if  these  values  change  from  tire  to  tire . 

l.  Deflection  cutoff  in  radii  (>-100  no  cutoff) 

This  is  a  new  parameter  not  used  in  the  current  CE  procedure, 
and  is  included  for  research  purposes  only.  If  the  deflection 
profile  under  a  wheel  load  is  plotted,  it  shows  that  the  theo¬ 
retical  deflection  only  gets  to  zero  at  a  very  long  distance 
from  the  center  of  the  tire.  In  reality,  this  deflection  pro¬ 
file  reaches  zero  at  considerable  shorter  offset  distances. 
Figure  7  graphically  illustrates  this  concept.  The  default 
value  of  100  is  used  for  routine  design  conditions. 

m.  Proceed  with  calculations?  (Y/N) 

Answer  yes  to  this  item  if  all  data  have  been  entered.  A  win¬ 
dow  showing  the  progress  of  the  calculations  is  displayed.  All 
input  parameters  must  be  entered  for  this  option  to  be  active. 
By  default,  the  program  will  always  create  an  output  data  file 
names  ESVLF.OUT.  The  results  in  this  file  may  be  saved  with 
another  name  using  one  of  the  options  from  the  main  menu. 

n.  Display  and  Print  Results?  (Y/N) 

After  the  calculations  have  been  performed,  the  user  can  dis¬ 
play  or  print  the  results.  If  no  calculations  have  been  per¬ 
formed,  a  window  opens  to  warn  the  user  to  perform  the 
calculations  first. 

o.  Back  to  main  menu?  (Y/N) 

This  option  takes  the  user  back  to  the  main  menu. 
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Y 


Y 


b.  With  X,  Y  symmetry 

Figure  5.  Example  grids  system  generated  by  the 
program  for  a  twin- tandem  gear 
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INPUT  PPRPMETERS  FOR  DRTR  FILE  :  [  EXAMPLE  .IN  ] 


SaM  Contact  Proaauro  for  all  tiros?  (V/N) . 

Saao  Radius  of  contact  arsa  for  all  tiros?  (V/N) 

Nuabor  of  tiros  in  goar  (108  Max) . 

Tiros 
Tiros 


Nuabor 
Dopth  I 
Grid  X- 
Grid  V- 
Syaaotr 
Syaaotr 


Dofloction  Cutoff  in  radii  (>=100  No  Cutoff) 

Procood  with  Calculations?  (V/N) . 

Display  and  Print  Results?  (V/N) . 

Back  to  Main  Menu?  (V/N)  . . 


> 


COORDINATES  OF 

TIRE  NUMBER 

1 

X-Coord.  (in) . 

. ► 

44.0 

V-Coord.  (in) . 

. ► 

0.0 

Contact  Prossuro  (psi) 

. ► 

200 

Radius  of  contact  aroa 

(sq  in)- • 

9.52 

Ok  ?  (V/N) . 

. ► 

V 

>  100 
N 
N 

>  N 


Figure  6.  Screen  for  the  input  of  tire  coordinates,  tire 
pressure,  and  radius  of  contact  area 


-  Actual  Deflection 

-  Theoretical  Deflection 

Figure  7 .  Theoretical  and  actual  deflection  profile 
under  a  circular  load 
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11.  Figures  8  and  9  present  the  example  problem.  Figure  10  shows  the 
screen  with  the  input  parameters  for  the  example.  The  arrow  keys,  page-up, 
and  page-down  keys  are  used  to  move  between  the  windows  input  fields  to  enter 
the  corresponding  values.  Once  all  parameters  are  entered  and  the  calcula¬ 
tions  performed,  the  option  to  display  or  print  the  results  can  be  selected. 
Figure  11  lists  the  results  in  a  tabular  form.  This  table  of  results  presents 
the  calculated  ESWL  factors  and  the  x  and  y  location  where  the  maximum  deflec¬ 
tion  occurs  for  a  specific  depth.  The  commands  available  are  displayed  at  the 
bottom  of  the  screen.  The  user  may  scroll  up  or  down  through  the  table  if 
necessary,  print  the  results,  or  get  back  to  the  data  entry  screen. 


Additional  Program  Features 


12.  The  program's  main  menu  contains  additional  features  to  save  the 
current  output  data  file  with  another  name,  view  previously  computed  results, 
and  load  previously  created  input  data  files. 

13.  Option  [2]  of  the  menu  allows  the  user  to  retrieve  an  input  data 
file  previously  created.  When  this  option  is  chosen,  a  list  of  the  available 
data  files  is  presented  from  which  the  user  may  select  one.  If  no  files 
exist,  the  program  warns  the  user  that  no  files  exist. 

14.  Option  [3]  allows  the  user  to  save  the  ESWLF.OUT  data  file  created 
from  a  previous  calculation.  A  window  then  asks  for  a  new  name  for  the  output 
data  file.  Every  time  calculations  are  performed,  the  program  automatically 
saves  the  results  to  the  ESWLF.OUT  file.  Because  this  operation  overwrites 
the  existing  ESWLF.OUT  file,  this  option  allows  the  user  to  save  results  from 
previous  calculations  to  a  file  with  another  name. 

15.  Option  [4]  allows  the  user  to  view  any  output  data  files  previously 
created  and  saved  with  option  [3].  A  list  of  available  output  data  files  is 
shown  on  the  screen  from  which  the  user  can  select  a  file  to  view. 
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Example  problem; 


Compute  the  equivalent  single-wheel  load  factors  (ESWLF)  for  one  side 
of  C-5A  main  ianding  gear.  Calculate  the  ESWL  down  to  100  in.  with 
increments  of  10  in.  Assume  no  deflection  cutoff  distance. 

Tire  Load  =  30,000  lb 

Tire  Contact  Area  =  285  sq.  in.  (measured) 

Tire  Radius  =  SQRT  (Tire  Contact  Area/3. 1415)  =  9.52  in. 


Y 


Figure  8.  Data  for  example  problem 


Solution: 


(1)  Tire  Coordinates 


re  No. 

X  (in) 

Y  (in) 

1 

-60.5 

0.0 

2 

-26.5 

0.0 

3 

26.5 

0.0 

4 

60.5 

0.0 

5 

-24.0 

65.0 

6 

24.0 

65.0 

7 

-60.5 

220.0 

8 

-26.5 

220.0 

9 

26.5 

220.0 

10 

60.5 

220.0 

11 

-24.0 

285.0 

12 

24.0 

285.0 

(2)  Number  of  depths  -  11  (10  increments  plus  surface  deflection) 

(3)  Depth  increment  -  10  in. 


(4)  Symmetry  along  Y-Axis 

(5)  No  Symmetry  along  X-Axis 

(6)  X  and  Y  grid  increments  -  2  in. 


Y-Axis 


Figure  9.  Tire  coordinates  and  grid  information  for  example  problem 
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INPUT  PARfiMETERS  FOR  ORTA  FILE  :  [  C5A-GEAR.IN  ] 


Saae  Contact  Pressure  for  all  tires?  (V/N) . ►  V 

Sane  Radius  of  contact  area  for  all  tires?  (V/N)->  V 

Number  of  tires  in  gear  (100  Hax) . ►  12 

Tire  Contact  Pressure  (psi) . ►  285.00 

Tire  Radius  of  contact  area  (in) . ►  9.52 

ESUL  Radius  of  contact  area  (in) . ►  9.52 

Number  of  computational  depths  (100  Hax) . ^  11 

Depth  Increment  (in) . ►  10.00 

Grid  X-Increment?  (in) . ►  2.00 

Grid  V-Increment?  (in) . ►  2.00 

Symmetry  along  X-Axis?  (V/N) . ►  N 

Symmetry  along  V-Axis?  (V/N) . ►  V 

Enter  or  change  Tire  Coordinates?  (V/N) . ►  N 

Deflection  Cutoff  in  radii  (>=100  No  Cutoff) . ►  100.00 

Proceed  with  Calculations?  (V/N) . ►  V 

Display  and  Print  Results?  (V/N) . ►  N 

Back  to  Main  Menu?  (V/N)  . ►  N 


a.  Input  parameters  for  example  problem 


Tire  H 

X-coord(in) 

V-coord( in) 

Radius( in) 

Press (psi) 

1 

-60.50 

.00 

9.52 

285.00 

2 

-26.50 

.00 

9.52 

285.00 

3 

26.50 

.00 

9.52 

285 . 00 

q 

60.50 

.00 

9.52 

285.00 

5 

-2*1.00 

65.00 

9.52 

285.00 

6 

2*1.00 

65.00 

9.52 

285.00 

7 

-60.50 

220.00 

9.52 

285.00 

8 

-26.50 

220 . 00 

9.52 

285.00 

9 

26.50 

220.00 

9.52 

285.00 

10 

60.50 

220.00 

9.52 

285.00 

11 

-2*1.00 

285.00 

9.52 

285.00 

12 

2*1.00 

285.00 

9.52 

285.00 

b.  C-5A  landing  gear  tire  coordinates 
Figure  10.  Example  problem  computer  screens 
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DEPTH(in) 

ESULF 

XMAX(in) 

VMftX(in) 

.08 

.891 

-26.50 

220.00 

18.00 

.152 

-26.58 

220.00 

20.00 

.197 

-28.58 

220.00 

30.00 

.2H7 

-28.58 

222.00 

140. 00 

.296 

-2‘4.5e 

22‘4.00 

50.00 

.3M*I 

-12.50 

226.00 

60.00 

.392 

-.50 

230.00 

70.00 

.‘433 

-.50 

230.00 

80.00 

.‘470 

-.50 

228.00 

90.00 

.503 

-.50 

228.00 

100.00 

.533 

-.50 

22‘4.00 

Figure  11 .  Table  showing  the  results  obtained 
from  the  example  problem 
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PART  IV:  CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 


16.  A  computer  program  for  the  calculation  of  the  equivalent  single¬ 
wheel  load  based  on  the  Boussinesq  solution  and  using  deflection  as  the  equiv¬ 
alency  parameter  was  presented.  The  following  conclusions  were  made: 

a.  Implementing  the  ESWL  concept  into  a  user  friendly  computer 
program  is  a  tedious  task  and,  although  the  concept  is  fairly 
simple,  the  numerical  analysis  for  its  solution  can  become  very 
complex. 

b.  The  computer  program  resulting  from  the  implementation  of  ESWLF 
concept  into  a  FORTRAN  code  produced  an  efficient  and  easy  to 
use  pavement  analysis  tool.  Sufficiently  accurate  results  for 
pavement  analysis  can  be  obtained  with  minimum  interaction 
between  the  user  and  the  program,  thus  hiding  the  user  from  the 
complex  mathematical  algorithms  involved  in  the  solution. 

c.  The  easy  to  use  interface  greatly  enhances  the  program  aesthet¬ 
ics  and  functionality,  and  eliminates  cumbersome  and  confusing 
data  entry  screens . 

Recommendations 

17.  Based  on  the  experience  gained  while  developing  and  using  this 
program,  it  is  recommended  that: 

a.  Grid  intervals  of  about  half  the  contact  area  radius  give  suf¬ 
ficiently  accurate  results. 

b.  Gear  symmetry  should  be  used  whenever  possible  to  reduce 
calculations . 
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APPENDIX  A:  COMPUTER  PROGRAM  SOURCE  CODE 


E;\£SWLF\ESULF.FOR  ■  Wed  Aug  21  10:18:08  1991 


SLARGE 

C  Equivalent  Single  Wheel  Load  Factor  program 
C  written  by:  Carlos  R.  Gonzalez 
C  Pavement  System  Division 

C  Geotechnical  Laboratory 

C  Waterways  Experiment  Station 

C  Vicksburg,  Mississippi 

(;«**»*».«.«.•.**«***•.**••••*•*••«••**•******••. 


INCLUDE  'FGRAPH.FP 
PROGRAM  ESWLFACTOR 

CHARACTER»12  FNAME,  OUTFILE,  WORKER,  FSPEC*40 
CHARACTER*3  EXT 
LOGICAL  QUIT 

CALL  C  F1LLSTR(FNAME, ■  •) 

CALL  C~FILLSTR<WORKER, '  ') 

CALL  cjiLLSTRCOUTFILE,  ■  ') 

CALL  c'nUHOFF 
CALL  Hello 
CALL  ScreenBk 
QUIT  =  .FALSE. 

DO  WHILEC. NOT. QUIT) 

CALL  MainMenul iPick) 

SELECT  CASE<iPick) 

CASEd) 

EXT  =  'IN' 

CALL  NewFNAME(WORKER,EXT) 

IF  (WORKER. N£.'  ')  THEN 
FNAME  =  WORKER 
CALL  ESWLF(FNAME) 

CALL  C  F I LLSTR (WORKER,'  ') 

END  I F 
CASE(2) 

FSPEC  =  '♦.IN' 

CALL  GetFile(WORKER, FSPEC) 

IF  (WORKER. NE.'  ')  THEN 
FNAME  =  WORKER 
CALL  C  FILLSTR(WORKER,'  ') 

END  IF 

IF  (FNAME. NE.'  ')  THEN 
CALL  ESWLF( FNAME) 

CALL  C  FILLSTR(WORKER,'  ') 

END  IF 
CASE(3) 

EXT  =  'OUT' 

CALL  NewFNAME(WORKER,EXT) 

IF  (WORKER. NE.'  ')  THEN 
CALL  SRESULTS(WORKER) 

CALL  C  FILLSTR(WORKER,'  ') 

END  IF 
CASE (4) 
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FSPEC  =  ••.OUT* 

CALL  GetFile(WORKER, FSPEC) 

IF  (WORKER. NE.'  *)  THEN 
OUT FILE  =  WORKER 
CALL  C  FILLSTR(WORKER.'  •) 
END  I F 

IF  (OUTFILE.NE.'  ')  THEN 
CALL  ShowRes(WORKER,OUTFILE) 
CALL  C  FILLSTR(WORKER,'  •) 
CALL  C~FILLSTR(OUTFILE,  ■  ') 
END  I F 
CASE(5) 

CALL  ShutDown 
QUIT  =  .TRUE. 

CASE  DEFAULT 
END  SELECT 
END  DO 

STOP  '  ' 

END 


SUBROUTINE  MainMenu( iPick) 


C 


CHARACTERV8  MA(9),  Title 


Hi  terns  =  9 

Title  =  '  ESWLF  PROGRAM  MAIN  MENU' 
MA  «  •  ' 


MACD  *  •  t1) 
MA(3)  =  •  C2] 
MA(5)  =  •  [3] 
MA(7)  =  •  [4] 
MA(9)  *  •  C5J 


Build  a  new  data  file' 

Load  an  existing  data  file  from  disk' 
Save  current  output  data  file  to  disk' 
View  an  output  data  file' 

Exit  program' 


CALL  C_CMENU(Nitems,MA,Title,iPick) 


RETURN 

END 


SUBROUTINE  NewFNAME<FNAME,EXT} 


CHARACTER*78  BA(1),  Title,  MSG 
CHARACTER  worker's,  FNAME*12,  EXT«3 
LCXilCAL  Escape,  exists,  VALID 
INTEGER  row,  col 


99  Ni  terns  =  1 

IF  (EXT. EQ. 'OUT')  THEN 
Title  «  'OUTPUT  DATA  FILE' 

BA(1)  *  '  Enter  name  of  output  data  file  t'.OUTIuuuul' 
ELSE 

Title  *  'NEW  DATA  FILE' 

BA(1)  B  ■  Enter  name  of  new  data  file  I'.INluOuui' 

END  IF 
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worker  =  '  ‘ 

CALL  ScreeriBk 

CALL  C_CBOX ( N i terns , BA , T i t I e, row, col , 8) 

77  CALL  C  GETSTRfrow, col, worker, 8, .TRUE. ..FALSE., .FALSE., 
&  ~  RR,1 I, lERR, Escape, lARROW) 

IF  (Escape)  THEN 
FRAME  =  ‘  ■ 

ELSE 

IF  (worker. EO.‘  ')  GOTO  77 
CALL  C  CHICFNAME(worker, VALID) 

IF  (.NOT. VALID)  THEN 
worker  =  '  ' 

GOTO  77 
END  IF 

FNAME  =  worker// '.'//EXT 
INQUIRE(FILE=FNAHE,EXIST=existS) 

IF  (exists)  THEN 
TITLE  =  'FILE  ALREADY  EXISTS!' 

MSG  =  '  Press  any  key  and  enter  a  new  file  name.' 
CALL  ERRMSG( TITLE, MSG) 

CALL  ScreenBk 
FNAME  =  '  ' 

GOTO  99 
END  IF 
END  IF 

CALL  ScreenBk 


RETURN 

END 

SUBROUTINE  GetFi te(FNAME, worker) 

CHARACTER  FNAME*12,  FSPEC*40,  worker*12 
CHARACTER*78  TITLE,  MSG 
c  LOGICAL  ERROR 

FSPEC  =  worker 
CALL  ScreenBk 

CALL  C_DIRFILES(FSPEC,NFILES, FNAME) 

CALL  ScreenBk 
CALL  C  CURSOROFF 
c  IF  (ERROR)  THEN 

c  TITLE  =  'DISK  OR  MEMORY  PROBLEMS’ 

C  MSG  =  '  CANNOT  READ  DIRECTORY! ' 

c  CALL  ERRMSG(TITLE,MSG) 

c  CALL  ScreenBk 

c  FNAME  =  '  ' 

c  ELSE 

IF  (NFILES.GT.O)  THEN 
•  Ok  to  return  FNAME 
ELSE 

TITLE  »  'FILES  NOT  FOUND' 

MSG  =  ’  NO  FILES  FOUND  MATCHING  SPEC:  '//FSPEC 
CALL  ERRMSG(TITLE,MSG) 

CALL  ScreenBk 
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FNAME  =  '  ' 
END  IF 

c  END  IF 

CALL  C_CURSORON 

RETURN 

END 


SUBRCXJTINE  SRESULTS< FNAME) 

CHARACTER* 12  FNAME 
CHARACTER*80  COUMMY 
LOGICAL  exists 

IDISK1  =  3 
I0ISK2  =  4 

INQUIRE(FILE='ESWLF.OUT*,EXIST=exists) 

IF  <.NOT.  exists)  THEN 
CALL  NoResults 
RETURN 
ELSE 

OPEN ( I D I SKI , F I LE= ' ESWL  F .OUT ' , STATUS= • OLD ' ) 
OPEN ( ID  I SK2 , F I LE=FNAME , STATUS= • NEW' ) 

DO  i=1,  1000 

REAO<I0ISK1,'(A)',EN0=1O)  COUMMY 
WRITE<IDISK2,'(A)')  CDUHMY 
END  00 

10  CLOSE(IOISKI) 

CLOSE<IOISK2) 

END  IF 

RETURN 


SUBROUTINE  ERRHSG(TITLE,MSG) 
CHARACTER*78  MSG,  BA(1),  Title 

Nitems  *  1 
BA(1)  s  MSG 
CALL  C  CURSOROFF 

CALL  C~CBOXFLASH(Nitenis, BA, Title) 
CALL  C"WAITF0RKEY 
CALL  C^CURSORON 

RETURN 

END 


SUBROUTINE  Hello 

CHARACTER  BA(11)*78  IBoxArray 
CHARACTER  Title*78 

BA  =  •  ■ 

Nlteois  *  11 

Title  ■  '  EQUIVALENT  SINGLE  WHEEL  LOAD  FACTOR  PROGRAM' 
BA<2)  =  '  written  by' 
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BA(4)  = 
BA(5)  = 
BA(6)  = 
BA{7)  = 
BAO)  = 
BA(9)  = 
BA(11)= 


Carlos  R.  Gonzalez* 

Pavement  System  Division* 
Geotechnical  Laboratory* 

US  Army  Engineer  Waterways  Experiment  Station* 
Vicksburg,  Mississippi* 

July  1991* 

Press  any  key  to  continue.* 


CALL  C_CURS0R0FF 
CALL  C'CLRSCR 

CALL  C  CBOX(Nitems,BA,Title,ir,ic,0) 
CALL  C'WAITFORKEY 
CALL  C'CURSORON 

RETURN 
END 
***** 

SUBROUTINE  Uait(MSG) 


CHARACTER  BA(1)*78  IBoxArray 
CHARACTER  Title* 78 
CHARACTER  MSG 


BA  =  *  ' 

N I  terns  =  1 
Title  =  MSG 

BA(1)  =  '  Please  wait...* 

CALL  C  CURSOROFF 

CALL  c“cB0X(Nitems,BA,Title,ir,ic,0) 

CALL  C_CURSORON 
CALL  ScreenBk 

RETURN 

END 

SUBROUTINE  ScreenBk 
CHARACTER  tp*80,  bt*80 

tp  =  'Equivalent  Single  Wheel  Load  FactorCESWLF)  Program* 
bt  =  *US  Army  Engineer  Waterways  Experiment  Station* 

CALL  C_SCREENBK(tp,bt,. FALSE.) 

RETURN 

END 

C******************************************************************** 

SUBROUTINE  ShutDown 


CALL  C_CLRSCR 

RETURN 

END 

(^*************************, 

SUBROUTINE  StatusBox 
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CHARACTER  BA(3)*78  'BoxArray 
CHARACTER  Title*78 

BA  =  '  ' 

N  Items  =  3 

Title  =  '  PERFORMING  CALCULATIONS’ 
CALL  C_CURSOROFF 

CALL  c“CBOXFLASH(Nitems, BA, Title) 
CALL  C_CURSORON 

RETURN 

END 


SUBROUTINE  STATUS( ICOUNT ,  I  TOTAL) 

INCLUDE  'FGRAPH.FO' 

RECORD  /rccoordy  s 
INTEGER*2  r,c 
INTEGER*2  duniny2,fg 
INTEGER*4  dummyA.bk 
CHARACTER  pct*4,  strng*80 

fg  =  GetTextColorO 
bk  s  GetBkColorO 

dumny2  =  SetTextColor(15)  Iwhite 

dumtrtyA  s  SetBkColor(4)  Iblue 

X  =  ICOUNT 

XT  s  I  TOTAL 

xpct  *  <X/XT)*100. 

ipct  =  xpct 

WRITE(pet,'(I3)')  ipct 
strng  =  pct//'X  COMPLETE' 
r  =  13 

c  =  <80-LEN_TRIM<strng))/2  ♦  1 
CALL  SetTextPosition<r,c,s) 

CALL  OutText(strng(1;LEN_TRIM(strng))) 
dainy2  =  SetTextColor(fg) 
dunmy4  =  SetBkColor(bk) 

return 

END 


SUBROUTINE  ESWLF(FNAME) 

INCLUDE  'FGRAPH.FD' 

CHARACTER  FNAME*12,  0UTFILE*12 
LOGICAL  QUIT,  CHANGES 
INTEGER  c 
INTEGER  jr,  jc 
INTEGER  row(20),  col 
LOGICAL  SAMEp,  SAMEr,  SAVE 

DIMENSION  X(100),  Y(100),  PRESS(IOO),  RAOIUSdOD),  Z(100) 
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CHARACTER*10  Field<20) 

DIMENSION  ICO0E(20) 

CHARACTER*20  Ch 
LOGICAL  Ok 

LOGICAL  Fexists 
CHARACTER  INTF'10 

LOGICAL  Escape 

CHARACTER  YESNO*1 ,  FNUM3*3,  FNUM7*7 
INTEGER»2  ITEM 

C  CALL  Uait( 'LOADING  PROGRAM') 

CALL  ReadFile(FNAME,NTIRES,X,Y, PRESS, RADIUS, RESW.NDEPTHS, 
&  DZ,GDX,GDY,1XSYM, I YSYM,RCUT, Fexists) 


SAMEp  =  .TRUE. 
SAMEr  =  .TRUE. 


NI terns  =  17 

IF  (.NOT. Fexists)  THEN 
DO  i  =  1,  NItetns 
DO  j=1,  10 
Field(i)(j:j)  =  •  • 

END  DO 
END  DO 
DO  i=1,  U 
ICOOE(i)  =  0 
END  DO 

NTIRES  =  1  !  set  some  defaults 

NDEPTHS  =  1 

PRESS(I)  =0.0 

RAOIUSd)  =  0.0 

RESW  =  0.0 

DZ  =  0.0 

GDX  =  1.0 

GOY  =  1.0 

RCUT  =  100.0 

Field(l)  =  'Y' 

Field(2)  =  'Y' 

Field(3)  =  '1' 

Field(«)  =  '0.0' 

Field(5)  =  '0.0' 

Fielded)  =  '0.0' 

Field(7)  =  '1' 

FieldCS)  =  '0.0' 

Field(9)  =  '1.0' 

FielddO)  =  '1.0' 

Fielddi)  =  'N' 

Fieldd2)  =  'N' 

Fieldd3)  =  'Y' 

Field(U)  =  '100' 

FielddS)  =  'N' 
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Field(16)  =  'N' 

Field(17)  s  'N' 

ICODE(I)  =  1 
IC00E(2)  =  1 
I  CODECS)  =  1 
1C00E(7)  =  1 

I  CODECS)  =  1 
IC0DE{9)  =  1 
ICODE(IO)  =  1 
ICODECll)  =  1 
1C0DE(12)  =  1 
ICODECIS)  =  1 
ICODE(14)  =  1 
CHANGES  =  .TRUE. 

ELSE 

DO  i=1,  NTIRES 

IF  (  PRESSCD.NE.PRESSCi)  )  SAMEp  =  .FALSE. 

IF  <  RAOIUSCD.NE.RADIUSCi)  )  SAMEp  =  .FALSE. 
END  DO 
DO  i=1,  14 
ICOOECi)  =  1 
END  DO 

CHANGES  =  .FALSE. 

DO  i=1,  NIteffls 
IF  (SAMEp)  THEN 
Field(l)  s  'Y' 

ELSE 

Fieldd)  *  'N* 

END  IF 

IF  (SAMEr)  THEN 
Field(2)  =  'Y' 

ELSE 

Field(2)  =  'N' 

END  IF 

WRtTE{INTF,'(nO)')  NTIRES 
READ<INTF,'<A10)')  Field(3) 

IF  <Field(1).EO.'Y')  THEN 
WRITE(INTF,'(F10,2)')  PRESSd) 
READ(INTF/(A10)')  Field(4) 

ELSE 

Field(4)  =  ' . . 

END  IF 

IF  <Field(2).EQ,'Y')  THEM 
WRITE<INTF,'(F10.2)')  RADIUSd) 
REA0(INTF/(A10)’)  Field(5) 

ELSE 

Field(5)  =  . 

END  IF 

WRITE(INTF,'<F10.2)')  RESW 
READ(INTF,'(A10)')  Field(6) 
WRITEdNTF,'(I10)')  NDEPTHS 
READ{INTF,'(A10)')  Field(7) 
WRITE<INTF,'<F10.2)')  OZ 
READ(INTF,'<A10)')  Field(8) 

IF  (DZ.EO.O.O)  Field(8)  =  'O-O’ 
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WRITEdNTF,  '(F10.2)')  GDX 
REA0{INTF,'{A10V)  Field(9) 

WRITEdNTF, '(FIO. 2)')  GDY 
READdNTF, '(A10)‘)  Field(IO) 

IF  dXSYM.EQ.1)  THEN 
Field(ll)  =  ‘Y* 

ELSE 

Field(ll)  =  'N' 

END  IF 

IF  (lYSYM.EQ.I)  THEN 
Fietd(12)  =  ‘Y* 

ELSE 

Fieldd2)  =  'N' 

END  IF 

FielddS)  =  'N' 

WRITEdNTF,  '(F10. 2)')  RCUT 
READdNTF,  ‘(AIO)')  Field(U) 

Field<15)  =  'N' 

Field(16)  =  'N' 

Fieldd?)  =  'N' 

END  DO 
END  IF 
DO  i=1,  17 

CALL  C  JUSTlFY(Fieldd),. FALSE.) 

END  DO 

CALL  DrawFormC jr, jc.FHAME) 
c  =  jc+1 
rowd)  =  jr 
DO  1=2,  NItems 

rowd)  =  rowd-1)  +  1 
END  DO 

DO  1=1,  N I  terns 

CALL  Wr1teFleld<row(i),c,Field(i)) 

END  DO 

ITEM  =  1 
QUIT  =  .FALSE. 

DO  WHILE  (.NOT. QUIT) 

SELECT  CASEdTEM) 

CASEd) 

99  c  =  jc+1 
col  =  Jc 
OK  =  .FALSE. 

DO  WHILE  (.NOT.  OK) 

YESNO  =  Fieldd)(1:1) 
chd:1)  =  YESNO 

CALL  C  GETSTR(row(1), col, YESMO.I,. true., .FALSE. ..false., 
£  RR,n,lERR, Escape,  (ARROW) 

SELECT  CASE  (YESNO(1:1)) 

CASE('y','Y') 

SAMEp  »  .TRUE. 
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OK  =  .TRUE. 

WRITE(IMTF,'(F10.2)')  PRESS(I) 

REA0(INTF,'(A10)')  FieldCA) 

CALL  C_JUSTIFY<Field(4),. FALSE.) 

CALL  UriteField(row(4),c,Field{4)) 

CASE('n','N') 

SAMEp  =  .FALSE. 

OK  =  .TRUE. 

Field(4)  =  . . 

CALL  WriteField(row(4),c,Field^4)) 

CASE  DEFAULT 
END  SELECT 
SELECT  CASE(IARROW) 

CASE(1:9) 

SELECT  CASE(YESNO) 

CASECY','y'.*N','n') 

OK  =  .TRUE. 

CASE  DEFAULT 
OK  =  .FALSE. 

END  SELECT 
CASE  DEFAULT 
OK  =  .FALSE. 

END  SELECT 
END  DO 

IF  (Ch(1:1).NE.YESN0)  THEN 
CHANGES  =  .TRUE, 

ICOOE(I)  «  1 
Fietd(1)<1:1)  »  yesNO 
ELSE 

Field<1)(1:1)  =  ch(1:1) 

END  IF 

CALL  C_ JUST lFY(Field(1),. FALSE.) 

CALL  WriteField(row<1),c,Field(1)) 

CASE (2) 

OK  *  .FALSE. 

DO  WHILE  (.NOT.  OK) 

YESNO  =  Field<2)(1:1) 
ch(1:1)  =  YESNO 

CALL  C_GETSTR( row(2) , col .YESNO, 1 , .TRUE. , . FALSE . , . FALSE . , 
'  RR, I l.IERR, Escape, lARROU) 

SELECT  CASE  (YESN0(1:1)) 

CASECy'.'Y') 

SAMEr  =  .TRUE. 

OK  =  .TRUE. 

WRITE(INTF,'(F10.2)')  RAOIUS(I) 

READ(INTF,'(A10)')  Field(5) 

CALL  C_JUSTIFY(Field(5),. FALSE.) 

CALL  WriteField<row<5),c,Field(5)) 

CASECn'.'N') 

SAMEP  »  .FALSE. 

OK  =  .TRUE. 

Field(5)  =  ' . . 

CALL  WriteField(row(S),c,Fieid'5T> 

CASE  DEFAULT 
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END  SELECT 
SELECT  CASE(IARROW) 

CASE(1:9) 

SELECT  CASE(YESNO) 

CASECY'/y'.'N'.'n') 

OK  =  .TRUE. 

CASE  DEFAULT 
OK  =  .FALSE. 

END  SELECT 
CASE  DEFAULT 
OK  =  .FALSE. 

END  SELECT 
END  DO 

IF  (ch(1:1).NE.YESNO)  THEN 
CHANGES  =  .TRUE. 

ICOD£(2>  =  1 
Field(2)(1;1)  =  YESNO 
ELSE 

Field(2)<1:1)  =  ch(1;1) 

END  IF 

CALL  C_JUSTlFY(Field(2), .FALSE.) 

CALL  UrjteFietd(row(2),c,Field(2)) 

CASE(3) 

OK  =  .FALSE. 

DO  WHILE  (.NOT.  OK) 

FNUM3  =  Field(3)<1:3) 

II  =  NTIRES 
IN  =  NTIRES 

CALL  C_GETSTR( row(3) ,col , FNUM3,3, . FALSE. , .TRUE . 
&  RR, 11, 1  ERR, Escape, I ARROW) 

IF  (lERR.EQ.O)  THEN 
IF  (II.GE.1  .AND.  ri.LE.100)  THEN 
OK  =  .TRUE. 

NTIRES  =  II 

Field(3)(1:10)  =  FNUM3 
IC00E<3)  =  1 
END  IF 
END  IF 

SELECT  CASE(IARROW) 

CASE(1:9) 

OK  =  .TRUE. 

CASE  DEFAULT 
OK  =  .FALSE. 

END  SELECT 
END  DO 

CALL  WrjteField(row(3),c,Field(3)) 

IF  (IN. NE. NTIRES)  THEN 
CHANGES  =  .TRUE. 

Field(13)(1:1)  =  'Y' 

END  IF 

CASE(A) 

IF  (SAMEp)  THEN 
DO  i=1, NTIRES 
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PRESS(i)  =  PRESS(I) 

END  00 
0(C  =  .FALSE. 

DO  WHILE  (.NOT.  OK) 

FNUM7  =  Fteld(4)<1:7) 

RR  =  PRESS(I) 

RN  =  RR 

CALL  C_GETSTR( row( A ) , col , FNUM7, 7, . FALSE ...  TRUE .,. TRUE 
S  RR, 11, lERR, Escape, lARROW) 

IF  (lERR.EQ.O)  THEN 

IF  (RR.GE.1.0  .AND.  RR.LE. 1000.0)  THEN 
OK  =  .TRUE. 

DO  j=1,  NTIRES 
PRESS(t)  =  RR 
END  DO 

Field{A)(1:10)  =  FNUM7 
ICOOE(A)  =  1 
END  IF 
END  IF 

SELECT  CASE (1 ARROW) 

CASE(1:9) 

OK  =  .TRUE. 

CASE  DEFAULT 
OK  =  .FALSE. 

END  SELECT 
END  DO 

CALL  WriteField(row(A),c,Field(A)) 

IF  (RN  .NE.  PRESSd))  CHANGES  =  .TRUE. 

END  IF 

CASE(5) 

IF  (SAHEp)  then 
DO  i=1, NTIRES 
RADIUS(f)  =  RADIUS(I) 

END  DO 
OK  =  .FALSE. 

DO  WHILE  (.NOT.  OK) 

FNUM7  =  Field(5)(1:7) 

RR  =  RADIUS(I) 

RN  s  RR 

CALL  C_GETSTR( row(5 ) ,col , FNUH7,7, . FALSE . , . TRUE . , .TRUE . 

RR, 1 1 , lERR, Escape, lARROW) 

IF  (lERR.EQ.O)  THEN 

IF  (RR.GE.0.1  .AND.  RR.LE. 1000.0)  THEN 
OK  =  .TRUE. 

DO  i=1,  NTIRES 
RADIUS(i)  <  RR 
END  DO 

Field{5)(1:10)  =  FNUM7 
ICOOE(S)  >  1 
END  IF 
END  IF 

SELECT  CASE (lARROW) 

CASE(1:9) 

OK  *  .TRUE. 
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CASE  DEFAULT 
OK  =  .FALSE. 

END  SELECT 
END  DO 
END  IF 

CALL  WriteFietd(row(5).c,Field<5)) 

CASE (6) 

OK  =  .FALSE. 

DO  WHILE  (.NOT.  OK) 

IF  (1COOE(6).EQ.O)  THEN 
FNUM7  =  Field<5)(1:7) 

RR  =  RAOIUS(I) 

ELSE 

FNUM7  =  Field(6)(1:7) 

RR  =  RESW 
END  IF 
RN  =  RR 

CALL  C_G£ TSTR( rou( 6 ) , co ( , FNUM7, 7, . FALSE ...  TRUE ., . TRUE . , 
&  RR, 1 1 , IERR,Escape, lARROW) 

IF  (lERR.EO.O)  THEN 

IF  (RR.GE.0.1  .AND.  RR.LE. 1000.0)  THEN 
OK  =  .TRUE. 

RESW  =  RR 

Field(6)(1:10)  =  FHUM7 
IC00E(6)  =  1 
END  IF 
END  IF 

SELECT  CASECIARROW) 

CASE(1:9) 

OK  =  .TRUE. 

CASE  DEFAULT 
OK  =  .FALSE. 

END  SELECT 
END  DO 

CALL  WriteField(row(6),c,Field(6)) 

IF  (RN  .NE.  RESW)  CHANGES  =  .TRUE. 

CASE (7) 

OK  =  .FALSE. 

DO  WHILE  (.NOT.  OK) 

FNUM3  =  Field(7)(1:3) 

II  =  NDEPTHS 
IN  s  II 

CALL  C_GETSTR(row(7),col ,FNUM3,3, .FALSE. , .TRUE. , . FALSE. , 
&  RR.II.IERR.Escape.IARROW) 

IF  (lERR.EQ.O)  THEN 
IF  (II.GE.1  .AND.  II.LE.100)  THEN 
OK  =  .TRUE. 

NDEPTHS  =  II 
Field(7){1:10)  =  FNUH3 
IC0DE(7)  =  1 
END  IF 
END  IF 

SELECT  CASE(IARROW) 
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CASE  (-1:9) 

OK  =  .TRUE. 

CASE  DEFAULT 
OK  =  .FALSE. 

END  SELECT 
END  DO 

IF  <IN  .NE.  NDEPTHS)  CHANGES  -  .TRUE. 

CALL  WriteField(row(7),c,F)eld(7)) 

CASE(8) 

OK  =  .FALSE, 

DO  WHILE  (.NOT.  OK) 

FNUM7  =  Field(8)(1;7) 

RR  =  02 
RN  =  RR 

CALL  C_GETSTfi( row(8) ,col , FNUM7, 7. . FALSE. , . TRUE . . . TRUE . , 
&  RR.n.IERR, Escape,  lARROW) 

IF  (lERR.EQ.O)  THEN 

IF  (RR.GE.0.0  .AND.  RR.LE. 1000.0)  THEN 
OK  =  .TRUE. 

02  =  RR 

IF  <D2  .EQ.  0.0)  THEN 
FNUM7  =  '0.0' 

Field(7)  =  M' 

NDEPTHS  =  1 
IC00E<7)  a  1 

CALL  WriteField(row(7),c,Field(7)) 

END  IF 

Field(8)(1:10)  *  FNUM7 
IC00E(8}  a  1 

end  if 

END  IF 

SELECT  CASE(IARROU) 

CASE(1:9) 

OK  a  .TRUE. 

CASE  DEFAULT 
OK  a  .FALSE. 

END  SELECT 
END  DO 
2(1)  a  0.0 
DO  ia2,  NDEPTHS 
Z{i)  a  2(i-1)  +  02 
END  DO 

CALL  Ur{teField(roH(8),c,Field(8)) 

IF  (RN  .NE.  02)  CHANGES  a  .TRUE. 

CASE(9) 

OK  a  .false. 

DO  WHILE  (.NOT.  OK) 

FNUH7  a  Field(9)(1;7) 

RR  a  GOX 
RN  a  RR 

CALL  C_GETSTR( row(9) , col , FNUH7, 7, . FALSE . , . TRUE . , .TRUE . , 
&  RR,II,IERR,Escape,IARROU) 

IF  (lERR.EO.O)  THEN 
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IF  (RR.GE.1.0  .AND.  RR .LE . 1000.0)  THEN 
0<  =  .TRUE. 

GDX  =  RR 

Field(9)(1:10)  =  FNUM7 
IC00E(9)  =  1 
END  IF 
END  IF 

SELECT  CASE(IARROW) 

CASE(1:9) 

OK  =  .TRUE. 

CASE  DEFAULT 
OK  =  .FALSE. 

END  SELECT 
END  DO 

CALL  UriteField(row(9),c,Field(9)) 

IF  (RN  .NE.  GDX)  CHANGES  =  .TRUE. 

CASEdO) 

OK  =  .FALSE. 

DO  WHILE  (.NOT.  OK) 

FNUM7  =  Field(10){1:7) 

RR  =  GDY 
RN  =  RR 

CALL  C_GET  STR ( row( 1 0 ) , co I , FHUM7, 7, . FALSE ...  TRUE ...  TRUE . , 
&  RR, 11, lERR, Escape, lARROW) 

IF  (lERR.EQ.O)  THEN 

IF  (RR.GE.1.0  .AND.  RR.LE. 1000.0)  THEN 
OK  =  .TRUE. 

GDY  =  RR 

Field(10)(l!l0)  =  FNUM7 
ICOOE(IO)  =  1 
END  IF 
END  IF 

SELECT  CASE(IARROW) 

CASE(1:9) 

OK  =  .TRUE. 

CASE  DEFAULT 
OK  =  .FALSE. 

END  SELECT 
END  DO 

CALL  WriteField(row(10),c,Field(10)) 

IF  (RN  .NE.  GDY)  CHANGES  =  .TRUE. 

CASEdI) 

OK  =  .FALSE. 

DO  WHILE  (.NOT.  OK) 

YESNO  =  Field(11)d;1) 
ch(1:1)  =  YESNO 

CALL  C_GETSTR ( rowd  1 ) , co I ,  YESNO,  1 , . TRUE . , .  FAL SE . , .  F AL SE . 
&  RR,II,lERR,Escape,IARROW) 

SELECT  CASE  (YESN0(1:1)) 

CASECy'.'Y') 

IXSYM  =  1 
OK  =  .TRUE. 

CASE( 'n' , 'N' ) 
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IXSYM  =  0 
OK  =  .TRUE. 

CASE  DEFAULT 

END  SELECT 

SELECT  CASE(IARROU) 

CASE(1;9) 

SELECT  CASE(YESNO) 

CAS£('Y',*y'.'N',*n') 

OK  =  .TRUE. 

CASE  DEFAULT 
OK  =  .FALSE. 

END  SELECT 
CASE  DEFAULT 
OK  =  .FALSE. 

END  SELECT 
END  DO 

IF  (ch(1:1).NE.YESNO)  THEN 
CHANGES  =  .TRUE. 

ICOOE(II)  =  1 
Field(ll)  =  YESNO 
ELSE 

Field(ll)  =  ch(1;1) 

END  IF 

CALL  C_JUST I F Y( F i e ld( 11 ) , . FALSE . ) 

CALL  UriteField(row(11),c,Field(11)) 

CASE(12) 

OK  *  .FALSE. 

DO  WHILE  (.NOT.  OK) 

YESNO  *  Field(12)(1:1) 
eh(1;1)  *  YESNO 

CALL  C_GETSTR(row(12),col,YESN0.1,. TRUE. ..FALSE. ..FALSE., 
‘  RR, II, lERR, Escape, lARROW) 

SELECT  CASE  (YESNO<1:1)) 

CASE('y'.'Y') 
lYSYM  =  1 
OK  *  .TRUE. 

CASECn'.'N') 
lYSYH  »  0 
OK  =  .TRUE. 

CASE  DEFAULT 
END  SELECT 
SELECT  CASE(IARROU) 

CASE(1:9) 

SELECT  CASE<YESHO) 

CASECY'.'y'.'N'.'n') 

OK  »  .TRUE. 

CASE  DEFAULT 
OK  «  .FALSE. 

END  SELECT 
CASE  DEFAULT 
OK  »  .FALSE. 

END  SELECT 
END  DO 

FIeld(12)  ■  YESNO 
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IF  (ch(1:1)  .NE.  YESNO)  THEN 
CHANGES  =  .TRUE. 

ICOOE(12)  =  1 
Field(12)  =  YESNO 
ELSE 

Field(12)  =  ch(1:1) 

END  IF 

CALI  C_JUSTlFY(Field( 12),. FALSE.) 

CALL  WriteField(row(12),c,Field{12)) 

CASE(13) 

OK  =  .FALSE. 

DO  WHILE  (.NOT.  OK) 

YESNO  =  Field(13)(1:1) 
ch(1:1)  =  YESNO 

CALL  C_GETSTR(row(13),COl,YESN0,1.. TRUE.,. FALSE.,. FALSE., 
&  RR, I 1,IERR, Escape, lARROW) 

SELECT  CASE  (YESNO(1:1)) 

CASECy',  'Y') 

SAVE  =  .TRUE. 

OK  =  .TRUE. 

CASE('n','N') 

SAVE  =  .FALSE. 

OK  »  .TRUE. 

CASE  DEFAULT 
END  SELECT 
SELECT  CASE(IARROW) 

CASE(1:9) 

SELECT  CASE(YESNO) 

CASECY'.'y'.'N'.'n') 

OK  s  .TRUE. 

CASE  DEFAULT 
OK  e  .FALSE. 

END  SELECT 
CASE  DEFAULT 
OK  r  .false. 

END  SELECT 
END  DO 

IF  {ch(1:1)  .NE.  YESNO)  THEN 
Field{13)  =  YESNO 
ELSE 

Field<13)  *  ch(1:1) 

END  IF 

CALL  C_JUSTIFY(Field(13),. FALSE.) 

CALL  WriteField(row(13),c,Field(13)) 

IF  (SAVE  .AND.  Field(13)(1:1)  .EO.'Y')  THEN 
IF  (Fexists)  THEN 
CHANGES  >  .FALSE. 

CALL  ShowTIres(NTIRES,X,Y, PRESS, RADIUS, 

&  SAMEp,SAMEr, Fexists, CHANGES) 

T  OVGES)  IC00E{13)  »  1 
Fie  .  3)  =  'N* 

CAL  v'r-;teField(row(13),c,Field(13)) 
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CALL  GetTi res(X,Y, PRESS, RADIUS, I.HTIRES, 

S  SAMEp,SAMEr,Fexists, CHANGES) 

Field(13)  =  'N' 

CALL  WriteField(row(13),c.Field(13)) 
IC00E<13)  =  1 
END  IF 

CALL  DrawForml jr, je,FNAME) 

END  IF 
C  =  jc  ♦  1 
DO  k=1,NItems 

CALL  WriteField(ro«(k),c,Field(k)) 

END  00 


CASE(U) 

OK  =  .FALSE. 

DO  WHILE  (.NOT.  OK) 

FNUM7  =  Field(14)(1:7) 

RR  =  RCUT 
RN  s  RR 

CALL  C_GETSTR ( row( 1 4 ) , co I , FNUM7, 7, . FALSE ... T RUE .,. TRUE . , 
i  RR.II.IERR, Escape, lARROW) 

IF  (lERR.EO.O)  THEN 
IF  (RR.GE.I  .AND.  RR.LE. 1000.0)  THEN 
OK  =  .TRUE. 

RCUT  =  RR 

Field<14)<1;10)  =  FNUH7 
IC00E(14)  »  1 
END  IF 
END  IF 

SELECT  CASE(IARROW) 

CASE(1:9) 

OK  =  .TRUE. 

CASE  DEFAULT 
OK  »  .FALSE. 

END  SELECT 
END  DO 

CALL  WriteField(row(14),c,Field(14)) 

IF  (RN  .NE.  RCUT)  CHANGES  =  .TRUE. 

CASE(IS) 

OK  =  .FALSE. 

SAVE  »  .FALSE. 

DO  WHILE  (.NOT.  OK) 

YESNO  =  'N' 

Ch(1:1)  »  YESNO 

CALL  C_GETSTR( row( 15 ) , co I .YESNO, 1 , . TRUE ...  FALSE ...  FALSE . , 
i  RR,II,IERR,Escape,IARROW) 

SELECT  CASE  (YESN0(1:1)) 

CASE('y','Y') 

SAVE  =  .TRUE. 

OK  »  .TRUE. 

DO  k»1,  14 

IF  (ICQOE(k}.EO.O)  THEN 
SAVE  ■  .FALSE. 
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YESNO  =  'N* 

EXIT 
END  IF 
END  00 

CASECn'.'N') 

SAVE  =  .FALSE. 

OK  =  .TRUE. 

CASE  DEFAULT 

END  SELECT 

SELECT  CASEdARROW) 

CASE {1:9) 

SELECT  CASE(YESNO) 

CASECY' ,  'y‘,  'N',  'n') 

OK  =  .TRUE. 

CASE  DEFAULT 
OK  =  .FALSE. 

END  SELECT 
CASE  DEFAULT 
OK  =  .FALSE. 

END  SELECT 
END  DO 

IF  (ch(1:1)  .NE.  YESNO)  THEN 
FielddS)  =  YESNO 
ELSE 

FielddS)  =  ch(1:1) 

END  IF 

CALL  C_JUSTIFY(Field(15),. FALSE.) 

CALL  WriteField(Pow(15),c, FielddS)) 

IF  (SAVE)  THEN 

CALL  SaveOata<FNAHE,NTIRES,X,Y,PRESS,RADIUS,RESW,NDEPTHS, 
4  OZ,GOX,GOY,IXSyM,dSYM,RCUT) 

CALL  GetESWLF(FNAME,NOEPTHS) 

CHANGES  z  .FALSE. 

CA  L  OrawForml jr, jc.FNAME) 
c  =  jc  ♦  1 
DO  k=1,NItenis 

CALL  Wr1teField(roM<k),e.Field(k)) 

END  DO 
END  IF 

CASEd6) 

OK  =  .FALSE. 

SAVE  =  .FALSE. 

Fieldd6)  =  'N' 

DO  WHILE  (.NOT.  OK) 

YESNO  =  Field(16)(1:1) 
ch(1:1)  =  YESNO 

CALL  C_GETSTR( roH( 16) , col , YESNO, 1 , .TRUE ...  FALSE ...  FALSE . , 
&  RR, II, lERR, Escape, lARROW) 

SELECT  CASE  (YESN0(1:1)) 

CASE('y','Y') 

IF  (CHANGES)  THEN 
SAVE  s  .FALSE. 

YESNO  =  'N' 
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CALL  ReCalc 

CALL  OrawFormCjr.jc, FRAME) 
c  =  jc  ♦  1 
DO  ks1,NItetns 

CALL  WriteField(row(k),c,Field(k)) 

END  DO 
ELSE 

SAVE  =  .TRUE. 

OK  =  .TRUE. 

END  IF 

CASECn'.'N') 

SAVE  =  .FALSE. 

OK  =  .TRUE. 

CASE  DEFAULT 

end  select 

SELECT  CASE(IARROU) 

CASE(1;9) 

SELECT  CASETYESNO) 

CASE('Y','y','N',*n') 

OK  =  .TRUE. 

CASE  DEFAULT 
OK  =  .FALSE. 

END  SELECT 
CASE  DEFAULT 
OK  =  .FALSE. 

END  SELECT 
END  DO 

IF  <ch(1;1)  .NE.  YESNO)  THEN 
Field(16)  »  YESNO 
ELSE 

Field(16)  =  ch(1:1) 

END  IF 

CALL  C.JUST I  FY( F i eld(16) , . FALSE . ) 

CALL  Wr i teF i e ld( row( 16) , c, F i eld( 16) ) 

IF  (SAVE)  THEN 
OUTFILE  s  'ESWLF.OUT' 

CALL  ShowRes(FNAHE, OUTFILE) 

CALL  OrawFormljr, je,FNAME) 
c  =  je  ♦  1 
00  k=1,NI terns 

CALL  WPiteField(row(k),c,Field(k)) 

END  DO 
END  IF 

CASE<17) 

OK  s  .FALSE. 

SAVE  *  .FALSE. 

Field(17)  a  'N' 

DO  WHILE  (.NOT.  OK) 

YESNO  a  Field(17)(1:1) 

Ch(1;1)  a  YESNO 

CALL  C_GETSTR( ro«( 17) , col , YESNO, 1 , . TRUE ...  FALSE ...  FALSE . , 
RR , 1 1 , 1 ERR .Escape, I ARROW) 

SELECT  CASE  (YESN0(1:1)) 
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CASECy'.'Y') 

SAVE  =  .TRUE. 

OK  =  .TRUE. 

CASECn'.'N') 

SAVE  =  .FALSE. 

OK  =  .TRUE. 

CASE  DEFAULT 

END  SELECT 

SELECT  CASE(IARROU) 

CASE<1:9) 

SELECT  CASEfYESNO) 

CASECY'.'y'.'N'/n') 

OK  =  .TRUE. 

CASE  DEFAULT 
OK  =  .FALSE. 

END  SELECT 
CASE  DEFAULT 
OK  =  .FALSE. 

END  SELECT 
END  DO 

IF  {ch(1:1)  .HE.  YESNO)  THEN 
Field(17)  =  YESNO 
ELSE 

Field(17)  =  ch(1:1) 

END  IF 

CALL  CJUST I FY ( F i eld<  1 7) , . FALSE .  ) 

CALL  UriteFi«(d(roM(17),c,Fie(d(17}) 

IF  (SAVE)  THEN 
QUIT  »  .TRUE. 

END  IF 

END  SELECT  !0F  ITEMS 
SELECT  CASE (I ARROW) 

CASE(A,8)  !up  &  left  arrow  keys 

IF(ITEH.EO.I)  THEN 
ITEM  =  17 
ELSE 

ITEM  =  ITEM  -  1 
END  IF 

CASE(2,5,6)  Idown  &  right  arrow  or  ENTER  keys 

IF  (ITEM. EQ. 17)  THEN 
ITEM  =  1 
ELSE 

ITEM  =  ITEM  ♦  1 
END  IF 

CASE(7,9)  'Home  or  PgUp  keys 

ITEM  =  1 
CASE{1,3) 

ITEM  =17  lEnd  or  PgDn  keys 

CASE(O)  'Esc  key 

(don't  move,  reset  field 

CASE  DEFAULT 
END  SELECT 
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END  DO  'CYCLE  UNTIL  EXIT  PROGRAM 


CALL  ScreenBk 


RETURN 

END 

(;**•****•*«**.•*•*.••*••••*••***••••** 

SUBROUTINE  DrawFonn(ir, ic.FNAKE) 

CHARACTER*78  BA(17),  Title 
CHARACTER  FNAME*12 


Title  =  'INPUT  PARAMETERS  FOR  DATA  FILE  :  [  '//FNAME//'  ]' 


BA  =  '  • 
NIteins  = 
BA(1)  =• 
BA(2)  =■ 
BA(3)  =■ 
BA(4)  =• 
BA(5)  =• 
BA(6)  =■ 
BA(7)  =• 
BA(8)  s' 
BA(9)  s' 
BA(10)=' 
BA<11)s' 
BA<12)*' 
BA(13)s' 
BA(U)s' 
BA<15)»' 
BA<16)s' 
BA(17)=' 


17 

Same  Contact  Pressure  for  alt  tires?  (Y/N)uuuuuuuiM' 
Same  Radius  of  contact  area  for  all  tires?  (Y/N)utM' 
Number  of  tires  in  gear  (100  Max)(AnjuuuuuuuuuuuuCnM' 
Tires  Contact  Pressure  (psi)uuuiUuuhjuuuuuuuOuuuuul' 
Tires  Radius  of  contact  area  (in)uuuuuuuuuuuuuuuuiM' 
ESUL  Radius  of  contact  area  (in)ihXiuuuuuuu(juuuuuutM' 
Nunber  of  computational  depths  (100  Hax)ijuuuuuuuu(A' 
Depth  Increment  (in)uuijuuOuiAnXnAAiuuuuuuuuuuuuuuu(M' 
Grid  X'lncrement?  (in)i^uuixjuuuuujutMuuiAiuOOuiX^ 

Symmetr/  along  X-Axis?  (Y/N)i^ijui>axAMXAjiMOuuuOuOuiM' 
Symmetry  along  Y-Axis?  (Y/N}0u0i>XxXUu0ijiUu0u0i>A^ 
Enter  or  change  Tire  Coordinates?  (Y/N)000u0uCi0000i' 
Deflection  Cutoff  in  radii  (>sioo  No  Cutoff )OOOOOOR' 
Proceed  with  Calculations?  (Y/N)i>jta30iMi>Uiji>j0u0(jij>' 
Display  and  Print  Results?  (Y/N)u0uuuu0uuuuuu0uuuul' 
Back  to  Main  Menu?  (Y/N)  i>>>juuOui>XjuuuuuuuuuuuuuuR ' 


CALL  C_CBOX(Nitems,BA,Title,irow,1col,10) 
ir  =  irow 
ic  s  ieol 


RETURN 

END 

Q*******«*******************«************************* 

SUBROUTINE  GetTires(x,y,PRESS,RADIUS,is,ie, 

&  SAMEp,SAMEr,FexistS, CHANGES) 

LOGICAL  SAMEp,  SAMEr,  Fexists,  CHANGES 

LOGICAL  Ok 
INTEGER  row(20),  COL 

CHARACTER  IHTF*3,  cn*3, INTF10*10 
CHARACTER  BA(5}*78 
CHARACTER  Title*78 


INTEGER  r.cjrje 
LOGICAL  SAVE 
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DIMENSION  X(100),  Y(IOO),  PRESS(IOO),  RADIUS(IOO) 
CHARACTER*10  Field(20) 

LOGICAL  QUIT,  Escape 
CHARACTER  FNUM7*7,  YESNO*1 
INTEGER*2  ITEM 

BA  =  ■  ' 

Nltems  =  5 

BA(1)  ='  X-Coord.  ( i n)uuuuuuOuuuOuuuuuuuuiM ' 

BA(2)  ='  Y-Coord.  < inTuuuuuuuuuuuuuuuuuuuui ' 

IF  (.NOT.SAMEp)  THEN 

BA(3)  ='  Contact  Pressure  (psi )uuuuuuuuuuui' 

END  IF 

IF  (.NOT.SAMEr)  THEN 

BA(4)  ='  Radius  of  contact  area  (sq  in)uuuV' 

end  if 

BA(5)  ='  Ok  ?  (Y/N)uuuuuuuuuuuuuuuuui' 

DO  i=is,  ie 
URITE(INTF,'<I3)')  i 
REA0(INTF,'(A3)')  cn 

Title  =• COORD I NATES  OF  TIRE  NUMBER  '//cn 
CALL  C_CBOX<Nitems,BA,Title,lr,ic,10) 
jr  =  ir 
jc  =  ic 
c  =  jc+1 
COL  s  jc 
row(l)  =  jr 
DO  k=2,5 

ro«(k)  s  row(k-l)  +  1 
END  DO 

IF  (.NOT.Fexists)  THEN 
Field  =  ‘  ‘ 

ELSE 

WRITE{INTF10,'(F10.2)')  X(i) 
REA0(INTF10,'(A10)')  Field(l) 

IF  (X(i)  .EQ,  0.0)  Field(l)  =  'O.O* 

CALL  C_JUSTIFY(Field<1),. FALSE.) 

CALL  WriteField<ro«(1),c,Field(1)) 
WRITE(INTF10,'(F10.2)')  Y(i) 
READ(INTF10,'(A10)')  Field(2) 

IF  (Y<i)  .EO.  0.0)  Field(2)  *  'O.O* 

CALL  C_ JUST  I F Y(  F  i e  I d( 2 ) , .  FALSE . ) 

CALL  WriteField(row(2),c,Ffeld(2)) 

IF  (.NOT.SAMEp)  THEN 
WRITE(INTF10,'(F10.2)')  PRESS(i) 
REAO(INTF10,'(A10)')  Field(3) 

CALL  C_JUSTIFY(Field(3),. FALSE.) 

CALL  WriteField(row(3),c,Field(3)) 

END  IF 

IF  (.NOT.SAMEr)  THEN 
WRITEdNTFlO, '(F10.2)')  RADIUS(i) 
READ(INTF10.'(A10)')  Field(4) 

CALL  C  JUSTIFY(Field(4), .FALSE.) 
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CALL  Uri teField<row(4),c.Field(4)) 

END  IF 

IF  (Fexists)  THEN 
Field(5)  =  'Y‘ 

ELSE 

Field(5)  =  'N* 

END  IF 
END  IF 

QUIT  =  .FALSE. 

IF  (Fexists)  THEN 
ITEM  =  5 
ELSE 
ITEM  =  1 
END  IF 

DO  WHILE  (.NOT.QUIT) 

SELECT  CASEtlTEM) 

CASE(I) 
r  a  jr 
c  =  jc+1 
OK  =  .FALSE. 

DO  WHILE  (.NOT.  OK) 

FNUM7  =  Field(1)(1:7) 

RR  =  X(i) 

CALL  C_GETSTR(row( 1 ) ,COL, FMUM7,7. . FALSE . , . TRUE . , . TRUE 
RR, 1 1 , lERR, Escape, lARROW) 

IF  (lERR.EQ.O)  THEN 

IF  ((RR.GE.-9999.)  .AMD.  (RR.LE.9999.))  THEM 
OK  a  .TRUE. 

X(i)  *  RR 

Field(1)(1;10)  a  FNUM7 

CALL  WriteField(ro«(1),c,Field(1)) 

CHANGES  a  .TRUE. 

END  IF 
END  IF 
END  00 

CASE(2) 

OK  a  .FALSE. 

DO  WHILE  (.HOT.  OK) 

FNUH7  a  Fie(d<2)<1:7) 

RR  a  Y(i) 

CALL  C_GETSTR( row(2), COL, FNUM7,7, . FALSE . , . TRUE . , . TRUE , 
RR, 1 1 , lERR, Escape, lARROW) 

IF  (lERR.EQ.O)  THEN 

IF  (<RR.GE.-9999.)  .AND.  (RR.LE.9999.))  THEN 
OK  a  .TRUE. 

Y(i)  a  RR 

Field(2)(1;10)  *  fnUM7 

CALL  WriteField(row(2>,c,Field(2)) 

CHANGES  a  .TRUE. 

END  IF 
END  IF 
END  DO 
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CASE(3) 

IF(.NOT.SAHEp)  THEN 
OK  =  .FALSE. 

00  WHILE  {.NOT.  OK) 

FNUH7  =  Field(3)(1;7) 

RR  =  PRESS(i) 

CAL L  C_GETSTR( row(3 ) , COL , FNUM7, 7, . FALSE . , . TRUE . , . TRUE . 
&  RR, 11, lERR, Escape, lARROW) 

IF  (lERR.EQ.O)  THEN 

IF  ((RR  .GE.  0.1)  .AND.  (RR.LE.1000. ))  THEN 
OK  =  .TRUE. 

PRESS{i)  =  RR 

Field(3)(1:10)  =  FNUM7 

CALL  WriteField{row(3),c,Field(3)) 

CHANGES  =  .TRUE. 

ENO  IF 
END  IF 
END  DO 
END  IF 

CASE (4) 

IF  (.NOT.SAMEr)  THEN 
OK  =  .FALSE. 

DO  WHILE  (.NOT.  OK) 

FNUM7  =  Field(4)(1:7) 

RR  =  RAOIUS(i) 

CALL  C_GETSTR( rOH(4 ) , COL , FNUM7, 7, . FALSE . , . TRUE . , . TRUE . , 
S  RR,II,IERR,Esc8pe,tARR0W) 

IF  (lERR.EQ.O)  THEN 

IF  ((RR.GE.  0.1)  .AND.  (RR.LE.1000.))  THEN 
OK  a  .TRUE. 

RAOIUSd)  a  RR 

Field{4)(1:10)  =  FNUM7 

CALL  WriteField(rOH(4),c,Field(4)) 

CHANGES  a  .TRUE. 

END  IF 
END  IF 
END  DO 
ENO  IF 

CASE (5) 

98  OK  a  .FALSE. 

Field(5)  a  <N' 

DO  WHILE  (.NOT.  OK) 

YESNO  a  Field(5)(1;1) 

CALL  C_GETSTR( row(5 ), COL , YESNO. 1 , .TRUE ...  FALSE ., . FALSE . 
S  RR , 1 1 , 1 ERR .Escape, I ARROW) 

SELECT  CASE  (YESN0(1:1)) 

CASECy'.'Y') 

SAVE  a  .TRUE. 

QUIT  a  .TRUE. 

OK  a  .true. 

CASECn'.'N') 

SAVE  a  .FALSE. 
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OK  =  .TRUE. 

END  SELECT 
END  DO 

CALL  WriteField{row(5),c,Fie(d(5)) 

END  SELECT 

SELECT  CASE(IARROW) 

CASE(8)  !up  arrow  key 

IF(ITEM.EQ.I)  THEN 
ITEM  =  5 
ELSE 

SELECT  CASEdTEM) 

CASE(5) 

IF  (SAMEr)  THEN 
ITEM  =  ITEM  -  1 
IF  (SAMEp)  ITEM  =  ITEM  -  1 
END  IF 
CASE (A) 

IF  (SAMEp)  ITEM  =  ITEM  -  1 
END  SELECT 
ITEM  =  ITEM  -  1 
END  IF 

CASE (2, 5}  Idown  arrow  or  ENTER  keys 

IF  (ITEM.EQ.S)  THEN 
ITEM  s  1 
ELSE 

SELECT  WSECITEM) 

CASE(2} 

IF  (SAMEp)  THEN 
ITEM  »  ITEM  ♦  1 
IF  (SAMEr)  ITEM  =  ITEM  ♦  1 
END  IF 
CASE(3) 

IF  (SAMEr)  ITEM  =  ITEM  ♦  1 
END  SELECT 
ITEM  »  ITEM  ♦  1 
END  IF 

^SE(7,9)  I  Home  or  PgUp  keys 

ITEM  »  1 
CASE(1,3) 

ITEM  »  5  !End  or  PgDn  keys 

CASE(O)  !Esc  key 

•don't  move,  reset  field 

CASE  DEFAULT 
END  SELECT 

END  DO  I  QUIT  LOOP 

END  DO  H  LOOP 

Fexists  =  .TRUE. 

RETURN 
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END 


SUBROUTINE  Wri teF i eld< i r, ic.strng) 


INCLUDE  'FGRAPH.FD* 

CHARACTER*(*)  Strng 
RECORD  /rccoord/  s 

INTEGER*2  duniny2,fg 
INTEGER*4  dunmyA.bk 

fg  =  GetTextColorO 

bk  =  GetBkColorO 

dutnmy2  =  SetTextColord  )  Iblue 

dumtnyA  =  SetBkColor(7)  Iwhite 


CALL  C_CURSOROFF 

CALL  SetTextPosition(ir,)c,s) 

CALL  OutText(strng) 

CALL  SetTextPosition(ir,ic,s) 
CALL  C  CURSORON 


dumny2  =  SetTextColor(fg) 
dummyA  =  SetBkColor(bk) 


RETURN 

END 


SUBROUT  I NE  ReadF  U  e( FNAME ,NT I  RES , X, Y . PRESS , RAD  I US , RESW, NDEPTHS , 
&  DZ,GDX,GOY,IXSYM,IYSYM,RCUT,FexistS) 


CHARACTER  FNAME*12 

DIMENSION  X(100),  Y(100),  PRESS(IOO),  RAD1USC100) 

CHARACTER  CDUMMY*80 
LOGICAL  Fexists 

IMQUrRE<FILE=FNAME,EXIST=Fexists) 

IF  {.NOT. Fexists)  RETURN 

CALL  C  OUNIT(IOISK) 

0PEN(1DISK,FILE=FNAHE) 

READdDISK.'CA)')  COUHMY 
READdOISK,*)  NTIRES 
IF  (NTIRES.  CT.  100)  THEN 

STOP  'ERROR:  Maximum  number  of  tires  should  be  <=  100.' 
END  IF 

READdDISK.'CA)')  COUMMY 
DO  i=1,  NTIRES 

READdOISK,*)  Xd),  Y(i),  PRESS(i),  RADIUS(i) 

END  DO 

READdOISK, '(A)')  CDUMMY 
READdOISK,*)  RESW 
READdOISK, '(A)')  COUMMY 
READdOISK,*)  NDEPTHS 
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REAOCIDISK, '(A)')  CDUMMY 
READdDISK,*)  OZ 
IF  (NOEPTHS  .GT.  100}  THEN 

STOP  '  ERROR:  Maximjn  nuniber  of  depths  should  be  <=  100.' 
END  IF 

READdDISK, '(A)'}  COUMHY 
READdDISK,*)  GOX,  GDY 
READdDISK, '(A)')  CDUMMY 
READdDISK,*)  IXSYM,  lYSYM 
READdDISK, '(A)')  CDUMMY 
READdDISK,*)  RCUT 

CLOSEdDISK) 

RETURN 

END 


SUBROUT 1 NE  Sa veOataf  FNAME , NT I RES, X, Y , PRESS , RAD I US , RESW, NDEPTHS , 
&  DZ,GDX,GOY, IXSYM, lYSYM, RCUT) 

CHARACTER  FNAME*12 

DIMENSION  XdOO),  Y(100),  PRESS(IOO),  RADIUSdOO) 

CHARACTER  S(20)*80 


S(1)='Number  of  tires  in  gear' 

S(2)='(X,Y)  tire  coordinates,  pressure,  radius  of  contact  area* 
S(3)='ESWL  Radius  of  contact  area' 

S(A)='Nut«ber  of  depths  (  if  *  1  ,  ESWLF  calculated  at  DZ)' 
S(5)='Depth  increment  (DZ)' 

S(6)*'DX,  OY  grid  increments' 

S(7)*'Gear  symmetry  along  X  *  Y  gear  axes  (0  if  not  symnetric)' 
S<8)*'0ef lection  Cutoff  Distance  in  radii  (if  >=  100  no  cutoff)' 

CALL  C  OUNITdDISK) 

OPENdOISK,FILE*FNAME) 

WRITEd0ISK,'{A)')  Sd) 

WRITEdDISK,'d3)')  NTIRES 
WRITEdOISK,'(A)')  S(2) 

DO  i>1,  NTIRES 

WRITEdOISK,'(4(F10.2))')  X(i),  Y(i),  PRESS(i),  RADIUSd) 

END  DO 

WRITEdDISK,'(A)')  S(3) 

WRITEd0ISK,'(F10.2)')  RESW 
WRITEd0ISK,'(A}'}  S(A} 

WRITEdOISK,'d3)')  NDEPTHS 
WRITEdD!SK,'(A)')  S(5) 

WRITEdDISK,'(F10.2}<}  OZ 
WRITEd0ISK,'(A)')  S(6) 

WRITEdOISK,'(2(F10.2)}'}  GOX.GOY 
WRITEdOlSK, '(A)')  S(7) 

WRITEdDISK,'(2l3)')  IXSYM,  lYSVM 
WRITEdOlSK, '(A)'}  S(8) 

WRITEdOlSK, '(FIO. 2)')  RCUT 

CLOSEdDISK) 
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RETURN 

END 

SUBROUTINE  ShowRes(FNAME,OUTFILE) 

CHARACTER*12  FNAME,  OUTFILE,  INFI  E 
LOGICAL  exists 
LOGICAL  OK,  GOPRINT 
CHARACTER*78  LA(IOO),  Title,  BA{5) 

CHARACTER*80  tp,  MSG,  COUMMY 

INQUIRE(FILE=OUTFILE,EXIST=exists) 

IF(. NOT. exists)  THEN 
CALL  NoResults 
RETURN 
END  IF 

CAL.  '  OUNIT(IDISK) 

OPENi.5iSK,FILE=OUTFILE) 

READdDISK, '(30X,A12)')  INFILE 
IF  (FNAME  .N£.  ‘  ■)  THEN 
IF  (INFILE. NE. FNAME)  THEN 
Ni terns  =  5 

Title  =  'CANNOT  PERFORM  REQUESTED  OPERATION!' 

BA  =  '  ' 

BA(1)  =  '  The  existing  output  data  file  does  not  match' 
BA(2)  s  '  the  data  in  the  input  data  file.  Re-calculate' 
BA(3)  =  '  the  ESULF  factors. ' 

BA(5)  =  '  Press  any  key  to  continue.' 

CALL  C_CURSOROFF 

CALL  C“CB0XFLASH(N  .ems, BA, Title) 

CALL  C_WAITFORKEY 
CALL  C~CURSORON 
CLOSE(TdISK) 

RETURN 
END  IF 
END  IF 

REAO(IOISK,'(A)')  Title 
READdDISK,  '(A)'  )  COUMMY 
DO  j=1,  100  'max  #  depths  »  100 
READdDISK, '(A)', END=10)  LA(j) 

END  DO 

10  Ni terns  =  j  -  1 
CLOSEdDISK) 

max  =  15  Imax  depths/window 
IF  (Nitems.LT.max)  max  =  Nitems 
Longest  =  1 
DO  i=1,  Nitems 

IF  (LEN  TRIH(LA(i)).GT. Longest)  Longest  =  LEN  TRIM(LA(i)) 
END  DO  “ 

IF  (LEN_TRIM(Title).GT. Longest)  Longest  =  LEN_TRlM(Title) 
icol  =  T+(80-Longest-2)/2 
irow  =  1+(25-max-4)/2 

TP  =  'ESWL  FACTOR  RESULTS' 
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MSG='[«l  PgUp  PgDn]  to  view.  {ENTER]  to  print,'// 

&  ‘  lEsc]  back  to  data  entry  screen* 

CALL  C  SCREENBK(TP,HSG.. FALSE.) 

OK  =  .FALSE. 

DO  WHILE(.NOT.OK) 

C  CALL  C_LIST<iroH,ieol,Nitems, max, LA, Title, iPick, Longest) 

CALL  C  SLIST<irow,icol,Nitenis, max, LA, Title, iPick) 
SELECT“CASE(iPiek) 

CASE(1:100)  (Print,  100=max  num.  depths 
OK  =  .TRUE. 

CALL  PrintRes(OUTFILE.GOPRINT) 

IF  (.NOT.GOPRINT)  OK  =  .FALSE. 

CASE(O)  (Esc 

OK  =  .TRUE. 

CASE  DEFAULT 
END  SELECT 

CALL  C  SCREENBKITP, MSG,. FALSE.) 

END  DO 

CALL  ScreenBk 


RETURN 

END 

r******************«**********M*****«*************' 

SUBROUT  I NE  ShowT i res( NT  I RES ,X , Y , PRESS , RAO 1 US , 

6  SAMEp,SAMEr,Fexists, CHANGES) 


C  INCLUDE  ' TOOLS. STR' 

CHARACTER  LA(100)*78,  Title»78 

CHARACTER*80  TP,  MSG 

CHARACTER*12  S1.  S2,  S3,  S4,  S5 

LOGICAL  SAMEp,  SAMEr,  Fexists,  CHANGES,  OK 

DIMENSION  X(100),  YdOO),  PRESS(IOO),  RAOIUS(IOO) 

Nitems  *  NTIRES 

C  1  12  24  36  48 

Title*  '  Tire#  X-coord(in)  Y-coord(in)'// 

&  '  Radius(in)  Press(psi)' 

DO  i=1,  Nitetns 
URITE{S1,'(I12)')  i 
WRITE(S2,'(F12,2)')  X(i) 

WRITE{S3,'(F12.2)')  Y<i) 

WRITE(S4,'(F12.2)')  RA0IUS<i) 

WRITE{S5,'{F12.2)')  PRESS(i) 

LA(i)  »  S1//S2//S3//S4//S5 
END  DO 

max  ■  IS  (max  depths/window 
IF  (Nitems. LT.awx)  max  *  Nitems 
Longest  *  1 
DO  i*1,  Nitems 

IF  (LEN  TRIM(LA(i)).GT. Longest)  Longest  *  LEN_TRIM(LA(i)) 
END  DO 

IF  (LEN_TRIM{Title).GT. Longest)  Longest  *  LEN_TRIM(Title) 
icol  ■  T+(80-Longest-2)/2 
irow  «  1+<25-max-4)/2 
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TP  =  'SELECT  TIRE  TO  CHANGE’ 

MSG='til  PgUp  PgDn]  to  mark,  [ENTER]  to  select,'// 

&  '  [Esc]  when  done' 

OK  =  .FALSE. 

DO  UHILE(.NOT.OK] 

00  i=1,  Nitems 

URITE(S1,'(I10)')  i 
WRITE(S2,'(F1Q.2)')  X(i) 

URITE(S3,'(F12.2)‘)  Y(i) 

WRITE(S4,'(F12.2)')  RADIUS(i) 

WRITEISS, '(F12.2)’)  PRESS(i) 

LA(i)  =  S1//S2//S3//S4//S5 
END  DO 

CALL  C_SCREENBK(TP, MSG,. FALSE.) 

C  CALL  C_LIST(iroH,icol, Nitems, max, LA, Title, iPick, Longest) 

CALL  C  SLIST(irow, icol, Nitems, max, LA, Title, iPick) 
SELECT“CASE(iPick) 

CASE(1:100)  ■100=max  num.  of  tires 

CALL  GetTires<X,Y, PRESS, RADIUS, iPick, iPick, 

&  SAHEp,SAHEr,Fexists, CHANGES) 

CASE(O)  !Esc 

OK  =  .TRUE. 

CASE  DEFAULT 
END  SELECT 
END  DO 

CALL  ScreenBk 

RETURN 

END 

SUBROUTINE  NoResultS 

CHARACTER  BA{6)*78,  Title»78 

Title  =  'OUTPUT  DATA  FILE  DOES  MOT  EXISTS!’ 

BA  =  '  ' 

Nitems  =  6 

BA(1)  z  '  No  output  data  file  has  been  generated.' 

6A(2)  =  '  Please,  make  sure  to  enter  your  data  and' 

BA(3)  s  '  calculate  the  ESWL  Factors  before  selecting' 

BA(4)  =  '  this  option. ‘ 

BA(6)  *  '  Press  any  key  to  continue.' 

CALL  C  CURSOROFF 

CALL  C~CBOXFLASH(Nitems,BA,Title) 

CALL  C'WAITFORKEY 
CALL  C^CURSORON 
CALL  ScreenBk 

RETURN 

END 

SUBROUTINE  ReCalc 
CHARACTER  BA(4)*78,  Title*78 
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BA  =  •  ■ 

Title  =  'PERFORM  CALCULATIONS  FIRST!' 

Nitefns  =  4 

BA(1)  =  '  Either  your  data  was  changed  or  no' 

BA(2)  =  '  ESWLF  calculations  have  been  performed.' 
BA(4)  =  '  Press  any  key  to  continue.' 

CALL  C  CURSOROFF 

CALL  C~CBOXFLASH(NitefflS. BA, Title) 

CALL  C'WAITFORKEY 
CALL  C“CURS0R0N 
CALL  ScreenBk 

RETURN 

END 

SUBROUT  I NE  Pr i ntResfOUT  F I LE .OkToPr i nt ) 

CHARACTER  OUTFILE*12,  INFILE*12 
LOGICAL  OK,  OkToPr int 
CHARACTER  0UMMY*80 

DIMENSION  X(100),  YdOO),  PRESS(IOO),  RADIUS(IOO) 
CHARACTER  BA(3)*78,  Title*78 

BA  =  '  ' 

NItems  =  3 

Title  »  'ABOUT  TO  PRINT  RESULTS' 

8A(1)  «  '  Please,  verify  the  printer  is  ready.' 
BA(3)  >  •  Press  (ENTER)  to  print,  (Esc)  to  CANCEL' 

CALL  C  CURSOROFF 
OK  s  .FALSE. 

OkToPrint  *  .FALSE. 

DO  WHILE  (.NOT .OK) 

CALL  C  CBOXFLASH(NiteiRS,BA,Title} 

CALL  c”gETASC(IKEY) 

SELECT  CASE(IKEY) 

CASE(13)  !  ENTER 

OK  *  .TRUE. 

OkToPrint  =  .TRUE. 

CASE(27}  !  Esc 

OK  =  .TRUE. 

OkToPrint  *  .FALSE. 

CASE  DEFAULT 
END  SELECT 
END  DO 

IF  (OkToPrint)  THEM 

C  CALL  C  OUNIT(LU) 

IDISK  =“3 

C  CALL  C  OUNIT(LU) 

IPRNT  *”4 

OPEN( IDISK, FILE«OUTFILE) 

READ(IDISK,'(30X,A12)')  INFILE 
CLOSE(IDISK) 

OPEN(IDISK,FILE>INFILE) 
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OPEN(IPRNT,FILE='LhT1 ■) 

WRITE(1PRNT,100)  iforn;  feed 
FORMATCIM 

WRITEdPRK'T, '(lx,*)')  'INPUT  PARAMETERS' 
READCIDISK, '(A)')  DUMMY 
WRITE(IPRNT,'(1X,A)')  DUMMY 
REAOdDISK, '(13)')  NTIRES 
WRITEdPRNT, '(1X.I3)')  NTIRES 
READCIDISK, '(A)')  DUMMY 
WRITEdPRNT, 'dX, A)')  DUMMY 
DO  i=1,  NTIRES 

READCIDISK, '(4(F10. 2))')  X(i),  Y(i), 

&  PRESS(i),  RADIUSd) 

WRITEdPRNT, '(1X.4(F10.2))')  X(i),  Y(i), 
&  PRESS(i),  RADIUSd) 

END  00 

READCIDISK, '(A)')  DUMMY 
WRITEdPRNT, '(IX, A)')  DUWY 
READCIDISK, '(F10. 2)')  RESW 
WRITEdPRNT, '(IX, Flo. 2)')  RESW 
READCIDISK, '(A)')  DUMMY 
WRITEdPRNT, '(IX, A)')  DUMMY 
READCIDISK, '(13)')  NDEPTHS 
WRITEdPRNT, '(IX, 13)')  NDEPTHS 
READCIDISK, '(A)')  DUMMY 
WRITEdPRNT, '(IX, A)')  DUMMY 
READCIDISK, '(F10. 2)')  02 
WRITEdPRNT, '<1X,F10.2)')  02 
READCIDISK, '(A)')  DUMMY 
WRITEdPRNT, '(lx, A)')  DUMMY 
READCIDISK, '(2(F10.2))')  a)X,GDY 
WRITEdPRNT,  ■(1X,2(F10.2))')  GDX,GDY 
READCIDISK, '(A)')  DUMMY 
WRITEdPRNT, '(IX, A)')  DUMMY 
READCIDISK, '(213)')  IXSYM,  lYSYM 
WRITEdPRNT, '(1X, 213)')  IXSYM,  lYSYM 
READCIDISK, '(A)')  DUMMY 
WRITEdPRNT, '(IX, A)')  DUMMY 
READCIDISK, '(A)')  DUMMY 
WRITEdPRNT, '(IX, A)')  DUMMY 
CLOSECIOISK) 

OPENdOISK,FILE=OUTFILE) 

WRITEdPRNT, '(IX, A)')  '  ' 

READCIDISK, '(A)')  DUMMY 
READCIDISK, '(A)')  DUMMY 
WRITEdPRNT, '(IX, A)')  DUMMY 
READCIDISK, '(A)')  DUMMY 
WRITEdPRNT, '(IX, A)')  DUMMY 
DO  i=1,  NDEPTHS 
READCIDISK, '(A)')  DUMMY 
WRITEdPRNT, '(IX, A)')  DUMMY 
END  DO 
CLOSECIOISK) 

CLOSECIPRNT) 

END  IF 
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CALL  C  CURSORON 

return" 

END 


SUBROUT  I NE  GetESWL  F ( FNAHE , NOEPTHS) 

implicit  REAL*8  (A-H,0-Z) 

CHARACTER  FNAME*12 

DIMENSION  X(IOO),  Y{100).  PRESS(IOO),  RADIUS(IOO),  2(100) 
DIMENSION  DEFMAXdOO).  XMAXD(IOO).  VMAXD(IOO) 

DIMENSION  ESUL(IOO),  ESWLF(IOO) 

CHARACTER  C0LIMMY*1 
LOGICAL  EXISTS 

CALL  C_OUNIT(LU) 
lOISK  =  LU 

C  CALL  C  OUNIT<LU) 

C  IFILE  =  LU 

I NQU I RE ( F I LE= FNAHE , EX I ST=EX 1 STS ) 

IF  (.NOT.  EXISTS)  THEN 
STOP  •  ERROR:  Data  file  does  not  exists.' 

END  IF 

IF  (NOEPTHS  .EQ.  1)  THEN 
OPEN(IFILE,FILE=<DEFL.OUT') 

END  IF 

0PEN(IDISK,F1LE»FNAME) 

REA0(1DISK,'(A)')  COUMMY 
REAOCtOIStC,*)  NTIRES 
IF  (NTIRES.  GT.  100)  THEN 

STOP  'ERROR:  Maximum  nuifcer  of  tires  should  be  <=  100.' 
End  if 

READ(IOISK,'(A)')  COUMMY 
DO  i=1,  NTIRES 

READdOISK,*)  X(i),  Y<i),  PRESS(i),  RADIUS(i) 

END  DO 

REAO(IDISK,'(A)')  COUMMY 
READdOISK,*)  RESW 
REAOdOISK,'(A)')  COUMMY 
READdOISK,*)  NOEPTHS 
READdOISK, '(A)')  COUMMY 
READdOISK,*)  02 
IF  (NOEPTHS  .GT,  100)  THEN 

STOP  '  ERROR:  Maximum  number  of  depths  should  be  <=  100.' 
END  IF 

IF  (02  .EO.  0.0)  NOEPTHS  =  1 
IF  (NOEPTHS. E0.1)  THEN 
2(1)  =  02 
END  IF 

00  i=2,  NOEPTHS 
2(i)  =  2(i-1)  ♦  02 
END  DO 

READdOISK, '(A)')  COUMMY 
READdOISK,*)  GOX,  GDY 
READdOISK, '(A)')  COUMMY 
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REAOdOISK,*)  IXSrW,  lYSVM 
REAOdOlSK, '(A)')  CDUHMY 
READdOISK,*)  RCUT 

CLOSEdOISK) 

IF  (GDX.EQ.0.0  .OR.  GOY.EO.O.O)  THEN 
STOP  '  ERROR:  Neither  GDX  or  GDY  can  be  eoual  to  zero  * 
END  IF  ^  ■ 

CALL  StatusBox 
CALL  C_CURSOROFF 

PI  =  3.  U1 592653589793238 
XMIN  =  999999.0 
VMIN  =  999999.0 
XMAX  =  -999999.0 
YMAX  =  -999999.0 
DO  1=1,  NTIRES 

IF  (X(i)  .LT.  XMIN)  XMIN  =  X{i) 

IF  (Yd)  .LT.  YMIN)  YMIN  =  Y(i) 

IF  (X(i)  .GT.  XMAX)  XMAX  =  X(i) 

IF  <Y(i)  .GT,  YMAX)  YMAX  =  Y(i) 

END  DO 

C  XMIN  =  XMIN  •  2.0*RA0IUS(1) 

C  YMIN  =  YMIN  -  2,0*RADIUS<1) 

C  XMAX  »  XMAX  +  2.0*RADIUS<1) 

C  YMAX  *  YMAX  ♦  2.0*RAOIUS<1 ) 

GXO  -  XMIN 

GYO  =  YMIN 

IF  (lYSYM  .EG.  1)  THEN 
XN  =  ((XMAX  •  XMIN)/2.0)/G0X 
ELSE 

XN  *  (XMAX  -  XMIN)/GDX 
END  IF 

IF  (IXSYM  .EG.  1)  THEN 
YN  =  ((YMAX  -  YMIN)/2.0)/GOY 
ELSE 

YN  =  (YMAX  -  YMIN)/GOY 
END  IF 

IF  (  XN  .EG.  INT(XN)  )  THEN 
NLX  =  XN  +  1.0 
ELSE 

NLX  =  XN  ♦  2.0 
END  IF 

IF  (  YN  .EG.  INT(YN)  )  THEN 
NLY  =  YN  +  1.0 
ELSE 

NLY  =  YN  +  2.0 
END  IF 

DO  i=1,  NDEPTHS 
OEFMAX(i)  *  -999999.0 
END  DO 
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ITOTAL  =  NDEPTHS*NLX*NLY*NTIRES 
1  COUNT  =  0 

DO  KZ=1,  NDEPTHS 
GY  =  GYO 
DO  KY=1,  NLY 
GX  =  GXO 
DO  XX=1,  NLX 
WZSUM  =0.0 
DO  KT=1,  HTIRES 

I  CUT  =  0  ! compute  defl.  at  cutoff  dist. 

IF  (RCUT.GE. 100.0)  ICUT  =  2  ! without  cutoff  distance 

88  CONTINUE 

IF  (ICUT.EO.O)  THEN 
RO  =  RCUT*RA0IUS(KT) 

ROFF  =  RO 
ELSE 

RO  =  SQRT(  (X(KT)-GX)**2  +  (Y«T)-GY)«*2  ) 

END  IF 

IF  (RO.LT. 0.000001)  R0=0.0 

R1  =  SORTC  Z(K2)**2  ♦  (R0-RAD1US{KT))**2  ) 

R2  =  SORT(  Z<KZ)**2  ♦  (R0*RADIUS(KT))**2  ) 

IF  (Z(K2)  .EQ.  0.0)  THEN 
IF  (RO.EO.O.O)  THEN 
WZ  =  (3.0*RA0IUS<KT)*PRESS(KT)/Z.0) 

GOTO  99 
END  IF 

IF  (RO.EQ.RAOIUS(KT))  THEN 
WZ  s  <3.0*RADIUS(ICT)*PRESS(ia)/PI) 

GOTO  99 
END  IF 

IF  (RO.GT.RAOIUS(KT))  THEN 
Sk  »  RADIUS<ICT>/RO 
skc  *  SQRT(1.0-sk**2) 
skc2  =  1.0-sk*»2 
CE  =  CEL(skc.1.00,1.O0.skc2) 

CIC  =  CEKskc,  1.00,1. 00,1. DO) 

XX  =  CE  -  (1.0-sk**2)*CK 

WZ  =  {(3.0*PRESS(ICT)*RAOIUS(ICT)/PI)*XX) 

GOTO  99 
END  IF 

IF  (RO.GT.0.0  .AND.  RO.LT.RADIUS(ICT))  THEN 
sk  =  RO/RADIUS<KT) 
skc  »  SORT(1.0-sk**2) 
skc2  =  1.0-sk**2 
CE  =  CEL(8kc,1.DO,1.DO,ske2) 

WZ  =  (3.0*RA0IUS<KT)«CE*PRESS{KT)/PI) 

END  IF 
ELSE 

IF  (RO  .EQ.  0.0)  THEN 

WZ  =  (3.0*PRESS<KT)*{RADIUS(ICT)*»2)/{2.0*R1)) 

ELSE 

COSV  =  {Z{KZ)**2  -  RAOIUS(l(T)**2  +  R0**2)/(R1*R2) 
skc  e  R1/R2 
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skc2  =  skc**2 

CE  =  CEKskc. 1.00.1. DO, skc2) 

CK  =  CEKskc, 1.00, 1.00,1. DO) 

XH  =  R2*CE 
XN  =  R1*CK 

U2  =  (3.0*PRESS(KT)*(XM-XM*COSV)/(2.0*Pn) 

END  IF 
ENO  IF 

99  continue 

IF  (ICUT.EQ.2)  THEN 
I  COUNT  =  I  COUNT  +  1 
CALL  STATUSdCOUNT.ITOTAL) 

ELSE 

IF  (ICUT.EQ.O)  THEN 
DOFF  =  W2 
ICUT  =  1 
GOTO  88 
ELSE 

IF  (RO.GE.ROFF)  THEN 
UZ  =  0,0 
ELSE 

UZ  =  W2  -  (DOFF*(RO/ROFF)) 

END  IF 

I  COUNT  =  I  COUNT  +  1 
CALL  STATUSdCOUNT.ITOTAL) 

ENO  IF 
END  IF 

WZSUH  UZSUM  +  UZ 
ENO  DO  j  KT 

C  IF  (NOEPTHS  .EQ.  1)  THEN 

C  URlTEdFlLE, '(3(F10,2))')  GX.GY, (-1  )*W2 

C  END  I F 

IF  (UZSUM  .GT.  OEFMAXCKZ))  THEN 
DEFMAXCKZ)  =  UZSUM 
XMAXD(KZ)  =  GX 
YMAXD(KZ)  =  GY 
END  IF 

GX  r  GX  +  GOX 
ENO  DO  I KX 

GY  =  GY  ♦  GOY 
END  00  ! KX 

ENO  DO  ! KZ 

C  IF  (NOEPTHS  .EO.  1)  CLOSEdFILE) 

GEARL  =0.0 
00  i=1,  NTIRES 

GEARL  =  GEARL  +  PRESS( i )*PI*(RAOIUS(i )**2) 

END  DO 

DO  i=1,  NOEPTHS 
R  =  SORT(  Z(i)**2  ♦  RESU»*2  ) 

ESUL(i)  =  2.0*R*PI»DEFMAX(i)  /  3.0 
ESULF(i)  =  ESUL(i)  /  GEARL 
ENO  DO 
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OPEN(IOISK,FILE='ESWLF.OUT*) 

WRITECIDISK.IOI)  FRAME 
WRITEdDlSK.lOO) 

DO  i=1,  NDEPTHS 

WRITECIDISK.IIO)  2(i).  ESWLF(i),  XHAXD(i),  YMAXD(i) 

C  URITEdOlSK.IIO)  Z(i),DEFMAXd>,ESWLF(i),XMAXD{i),YMAXD{j) 

END  DO 
CLOSEdDISK) 

101  FORMAT (IX, 'RESULTS  FOR  INPUT  DATA  FILE:  •.A12) 

C  +  ♦  +  +  + 

C  1 2345678901 2345678901234567890123456789012345678901234567890 

100  F0RMAT(2X,'  DEPTHdn)  ESWLF  XMAXdn)  YMAX(in)', 

^  j  2x  - ...... - - - .1^ 

110  FORMAT(Fl6.2,2X,F10.3,2X,2(F10.2)) 

C  100  FORMAT<2X,'  DEPTH(in)  WZ*Etn(tn*psi )  ESWLF', 

C  S'  XMAX(in)  YMAX(in)', 

C  S  . . . 

C  S' . ■) 

C  110  FORMAT(2(F10.2,2X),F10.3,2X,2{F10.2)) 

CALL  C  CURSORON 

return" 

END 


Complete  elliptic  integral  solution  (1st  and  2nd  kind)  from: 
'NUMERICAL  REClPeS(FORTRAN)* ,  CAMBRIDGE,  REPRINT  1990 
pp.  187-188 

C******************************************************************** 


FUNCTION  CcL(OQC,PP,AA,BB) 

IMPLICIT  REAL*8  (A-H,0-Z) 

PARAMETER  (CA=. 000001,  PI02=1 .5707963268) 

IF  (OQC  .EQ.  0.0)  THEN 
STOP  '  Failure  in  CEL  furwtion.' 

END  IF 

QC  >  ABS(OOC) 

A  s  AA 
B  r  BB 
P  =  PP 
E  =  QC 
EM  =  1.0 

IF  (P  .GT.  0.0)  THEN 
P  =  SQRT(P) 

B  =  B/P 
ELSE 

F  =  QC*QC 
Q  =  1.0  -  F 
G  =  1.0  -  P 
F  =  F  -  P 
0  =  O^tB-A^P) 
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P  =  SQRT<F/G) 

A  =  (A-B)  /  G 
B  =  -Q/(G*G*P)  +  A»P 
END  IF 
1  F  =  A 

A  =  A  +  B/P 
G  =  E/P 
B  =  B  ♦  F*G 
B  =  B+B 
P  =  G+P 
G  =  EM 
EM  =  OC+EM 

IF  (A8S(G-0C)  .GT.  G»CA)  THEN 
QC  =  SQRT(E) 

QC  =  OC+OC 
E  =  QC*EM 
GOTO  1 
END  IF 

CEL  =  PI02*(B+A*EM)/(EM*(EM+P)) 

RETURN 

END 
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